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Abstract 

Ultracold atoms uniformly filling an optical lattice can be treated like an artificial crystal. 
An implementation including the atomic occupation of a single excited atomic state can 
be represented by a two-component Bose-Hubbard model. Its phase diagram at zero tem- 
perature exhibits a quantum phase transition from a superfluid to a Mott insulator phase. 
The dynamics of electronic excitations governed by electrostatic dipole-dipole interactions 
in the ordered regime can be well described by wave-like collective excitations called ex- 
citons. Here we present an extensive study of such excitons for a wide range of optical 
lattice geometries and of different dimensionality including boundary effects in finite lat- 
tices. As they are coupled to the free space radiation field, their decay depends on the lattice 
geometry, polarization and lattice constant. Their lifetimes can vary over many orders of 
magnitude from metastable propagation to superradiant decay. Particularly strong effects 
occur in one dimensional atomic chains coupled to tapered optical fibers where free space 
emission can be completely suppressed and only a superradiant interaction with the fiber 
modes takes place. We show that coherent transfer of excitons among spatially separated 
optical lattices can be controlled and represents a promising candidate for quantum infor- 
mation transfer. For an optical lattice within a cavity the excitons are coupled to cavity 
photons and the resulting collective cavity QED model can be efficiently formulated in 
terms of polaritons as elementary excitations. Their properties are explicitly calculated for 
different lattices and they constitute a non-destructive monitoring tool for important sys- 
tem properties as e.g. the atomic quantum statistics. Even the formation of bound states 
and molecules in optical lattices manifests itself in modified polariton properties as e.g. an 
anisotropic optical spectrum. Partial dissipation of the exciton energy in the lattice leads to 
heating, which can be microscopically understood through a mechanism transferring atoms 
into higher Bloch bands via a resonant excitation transfer among neighboring lattice sites. 
The presence of lattice defects like vacancies in the Mott insulator induces a characteristic 
scattering of polaritons, which can be optically observed to monitor the lattice integrity. Our 
models can be applied to simulate and understand corresponding collective phenomena in 
solid crystals, where many effects are often masked by noise and disorder. 
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1 Introduction 

Trapping and manipulating ultracold neutral atoms in an array of optical poten- 
tials opened a vast new research area in the field of quantum fluids and solid state 
modeling (Bloch et al., 2008). Demonstrating the superfluid to Mott insulator quan- 
tum phase transition of interacting bosons in an optical lattice laid the foundation 
to study and simulate strongly correlated quantum many-body systems with a big 
impact on condensed matter physics (Lewenstein et al., 2007). Ultracold atoms in 
optical lattices are now one of the most flourishing directions of experimental and 
theoretical quantum physics (Jaksch et al., 1998; Greiner et al., 2002; Morsch and 
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Oberthaler, 2006). Besides fundamental physics issues on entanglement, measure- 
ment and decoherence (Mekhov and Ritsch, 2012), more and more applications 
appear for the quantum simulations for a wide range of puzzling effects in solid 
state physics. 

A BEC of a dilute bosons, in which a very large number of atoms occupy the same 
quantum state of motion, was experimentally realized by several groups (Dalfovo 
et al., 1999). Optical control of the quantum state of such BEC opened the door to a 
new generation of experiments (Metcalf and van der Straten, 1999). In particular the 
loading of the BEC into a standing wave formed by pairs of counter propagating far 
detuned laser beams creates a synthetic controllable periodic lattice (Greiner et al., 
2002; Spielman et al., 2007). Trapped atoms still can tunnel from site to site, while 
interacting and repelling each other when occupying the same site. The stronger 
the laser field, the deeper the lattice and the slower the hopping rate of the atoms, 
while at the same time the on-site interaction between the atoms becomes stronger. 
Thus the laser light allows to switch between a weakly interacting gas, where atoms 
form a superfluid, to a strongly correlated regime, where atoms avoid each other and 
localize into individual sites forming a perfectly regular array called Mott insulator 
phase. As this crossover happens at T = as a function of an external parameter, 
it is called a quantum phase transition (Fisher et al., 1989; Sheshadri et al., 1993; 
Sachdev, 1999). It is theoretically well described by a Bose-Hubbard model (Jaksch 
et al., 1998) and soon was experimentally demonstrated (Greiner et al., 2002). An 
increasing degree of control of experimental parameters in optical lattices has been 
demonstrated up to the level of addressing and controlling single atoms at a specific 
lattice site (Sherson et al., 2010; Bakr et al., 2010; Weitenberg et al., 201 1). 

In a parallel development the dynamics of the electromagnetic field modes in a 
high-Q optical cavity resonantly coupled to quantized excitations of atoms has been 
extensively studied. Many phenomena discussed as Gedanken-experiments to ex- 
hibit the quantum physics of coupled systems and measurements were actually per- 
formed experimentally and lead to the recent 2012 physics Nobel prize (Haroche 
and Raimond, 2006; Walther et al., 2006). The first experimental realizations of 
cavity QED used single Rydberg atoms in superconducting microwave resonators. 
Later ground state alkali atoms in high-Q optical resonators allowed to reach the 
strong coupling limit as well, where light matter coupling dominates decoherence 
(Ye et al., 1999; Pinkse et al., 2000). Other implementations using solid state tech- 
nology followed more recently (Kavokin and Malpuech, 2003). 

It is of course a natural development to combine these two threads and investi- 
gate ultracold atomic lattices within optical resonators (Maschler and Ritsch, 2005; 
Zoubi and Ritsch, 2007; Mekhov et al., 2007; Larson et al., 2008; Ritter et al., 
2009). As a major experimental step in this direction, several groups (Brennecke 
et al., 2007; Colombe et al., 2007; Slama et al., 2007) achieved operations of a BEC 
of trapped ultracold atoms within an optical high-Q cavity deeply in the strong res- 
onant coupling regime. It seems thus that no major technical obstacles prevent the 
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implementation of an optical lattice within an optical cavity. An alternative route 
for realizing an optical lattice within a cavity is based on a co-planar waveguide 
microwave resonator resonantly coupled to the collective hyperfine state of a BEC 
of ultracold atoms in a wire surface trap fabricated on a superconducting atom chip 
(Folman et al., 2002). This system is predicted to allow the experimental study 
of magnetic dipole interactions in the strong coupling regime (Verdu et al., 2009; 
Imamoglu, 2009). 

Placing ultracold atoms within optical cavities leads to a wealth of interesting dy- 
namical phenomena already in the dispersive regime, where the light frequency is 
far detuned from the atomic resonance. This involves a dynamic interplay of light 
forces, non-local couplings and field dynamics (Ritsch et al., 2012). In this review 
instead we concentrate on resonant interaction of quantized light modes and trapped 
atoms in the Mott phase, which we treat as a new type of matter closely related to an 
atomic crystal with widely controllable lattice geometry (Zoubi and Ritsch, 2007). 
While in some general cases we consider coupling of the light to motional degrees 
of freedom as well, our main emphasis is on collective internal excitations dynam- 
ics. This reduced description already allows to study many solid state effects in a 
very generic form. The wide and precise controllability of different lattices and cou- 
pling parameters allows to change the number of atoms per site, the geometry, and 
the symmetry of the system. These facts give the present system a decisive advan- 
tage over conventional solid state crystals where parameters are hardly tunable and 
particular effects cannot be easily singled out. Such toy implementations provide to 
us a more detailed and deeper understanding of solid state effects. Mathematically 
one extends the common Bose-Hubbard atom model to include atoms in excited 
states formally described as a second kind of bosons within the optical lattice. The 
system is then represented by a two-component Bose-Hubbard model exhibiting a 
much richer phase diagram (Zoubi and Ritsch, 2009d). 

Optical lattices allow us to go beyond the weakly interacting regime of Bose gas 
into strongly correlated systems of solid state physics (Bloch et al., 2008; Lewen- 
stein et al., 2007). Including the internal atomic level structure bears a strong anal- 
ogy to the excitonic dynamics of molecular crystals, described by "Frenkel exci- 
tons" (Agranovich, 2009; Davydov, 1971). The electronic excitations delocalize in 
the lattice due to electrostatic interactions, e.g. resonant dipole-dipole interactions. 
By exploiting the lattice symmetry one is lead to introduce wave-like collective 
electronic excitations called "excitons" as fundamental system excitations (Zoubi 
and LaRocca, 2005a). For this expansion to make sense, the transfer time has to 
be faster than the exciton lifetime, so that collective states of electronic excitations 
dominate the electrical and optical system properties. Excitons thus provide an ideal 
basis to describe resonant effects in optical lattices. In fact the atomic recoil heating 
due to emission and absorption of a photon is significantly suppressed as the mo- 
mentum is shared among all atoms in the lattice and strongly suppress diffusion and 
heating in the Mott insulator. We present studies of the properties and dynamics of 
excitons for several generic lattice types of different geometries and dimensional- 
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ity, starting from finite chains, infinite chains to stacks of two dimensional planar 
lattices (Zoubi and Ritsch, 2007, 2009c, 2010b, 2011a,b). Boundary conditions of 
finite systems give rise to standing wave excitons, while extended systems produce 
propagating wave solutions. As a special case a Mott insulator with two atoms per 
site can also possess dark on-site localized excitons decoupled from the remaining 
lattice sites (Zoubi and Ritsch, 2008b). 

The energy derealization of excitons among many particles amounts to the occur- 
rence of intrinsic entanglement. To generate and preserve significant entanglement 
requires sufficient isolation from the environment in order to bring decoherence and 
dissipation to a minimum, which is a big challenge for optically excited systems. 
Here spontaneous emission due to the coupling to free space radiation is a ma- 
jor source of dissipation and decoherence (Loudon, 2000). Excitons in multi-atom 
systems exhibit significantly modified spontaneous emission compared to a sin- 
gle atom (Ficek and Tanas, 2002), which we have extensively studied for different 
geometries and dimensions (Zoubi and Ritsch, 2010b, 201 la; Zoubi, 2012a; Oster- 
mann et al., 2012). While often excitons exhibit enhanced coupling to excitation 
lasers and vacuum modes and fast or superradiant decay, in many cases one can get 
long lived or even dark and metastable excitations, which can only decay optically 
into a narrow set of wave-vectors analogous to Bragg angles. If the lattice spacing 
is smaller than the transition wavelength, the exciton wave vector can be larger than 
the free space photon wave-vector, so that no optical decay is possible. Neverthe- 
less one can still engineer coherent transfer of such excitons among different lattice 
regions or different optical lattice planes (Zoubi and Ritsch, 2010c, 201 lb, 2012b). 
This effect should have important applications in optical information storage and 
transfer based on large ensembles. Interestingly, while the individual dipole-dipole 
energy transfer has an inverse cubic dependence on the distance, collective states 
within one lattice plane only couple with an exponentially decreasing strength to 
neighboring planes, which generally allows simplification to a nearest neighbor 
coupling model. 

Tailored light-matter interaction is of great importance to fundamental physics tests 
as well as high end technical applications (Cohen-Tannoudji et al., 1992; Metcalf 
and van der Straten, 1999; Meystre, 2001) and in particular for the development 
of quantum information technologies (Pinkse et al., 2000; Ye et al., 1999). New 
frontiers of cavity QED involving ultracold atoms trapped in intra-cavity optical 
lattices form a natural, but only very recently opened, research direction (Zoubi 
and Ritsch, 2007, 2008b, 2009c). As above, we concentrate here on atoms with an 
internal electronic transition close to resonance with a cavity mode, while the opti- 
cal lattice is generated by external classical fields which are off-resonance to any of 
the internal atomic transitions forming the dynamical part of the system. The case 
of cavity generated optical potentials also exhibits a wealth of interesting phenom- 
ena (Ritsch et al., 2012). Here we generalize the two-component Bose-Hubbard 
model to include cavity photons and their coupling to electronic excitations of the 
cold atoms (Zoubi and Ritsch, 2009d). The phase diagram including superfluid and 
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Mott regions is affected by the coupling to cavity photons which mediate long range 
interactions and enhance tunneling via field fluctuations. Even though using a rather 
simple model we see that the quantum phase transition to a perfectly ordered Mott 
state occurs at much deeper optical lattices than in free space. 

In cavities, excitons and photons are coherently coupled to form new quasi-particles 
called polaritons (Agranovich, 2009; Kavokin and Malpuech, 2003; Zoubi and 
LaRocca, 2005a). As the cavity mirrors exhibit a finite transition probability, the 
detection of the emitted light provides us with detailed information about the sys- 
tem properties (Zoubi and Ritsch, 2007). The corresponding optical spectra give a 
largely nondestructive observation tool for different quantum phases with minimal 
back-action (Mekhov and Ritsch, 2012). As an additional example the formation of 
molecules in the optical lattice introduces optical anisotropy in our artificial crystal. 
This manifests itself in the spectra through the cavity photon polarization mixing 
(Zoubi and Ritsch, 2009b). Again our studies range from finite lattices inside op- 
tical cavities of spherical mirrors up to two dimensional lattices between planar 
cavity mirrors. As an interesting property for the search of lower temperatures, the 
elastic scattering of excitons and polaritons by the defects or vacancies in the Mott 
insulator phase can be observed via cavity transmission spectra without destroying 
the lattice state as it occurs in high aperture microscope detection via fluorescence 
(Zoubi and Ritsch, 2008a). Similarly atoms excited to higher Bloch bands change 
the dynamics of excitons and polaritons (Zoubi and Ritsch, 2009a, 2010a). 

Optical properties of a gas of ultracold atoms with coherently dressed two-level and 
three-level systems in a Mott insulator phase of an optical lattice have been stud- 
ied by Bariani and Carusotto (2008); Carusotto et al. (2008). Cooperative effects 
of ultracold molecules in optical lattices are also treated (Kuznetsova et al., 2012). 
Recently Rydberg atoms become of large interest for optical lattices. Rydberg ex- 
citations in Bose-Einstein condensates of Rubidium atoms loaded into quasi-one- 
dimensional traps and in optical lattices experimentally have been realized (Viteau 
et al., 201 1 ; Anderson et al., 201 1). Moreover, the formation of excitons in a regular 
and flexible chain of Rydberg atoms as introduced through the interplay between 
excitonic and atomic motion was proposed theoretically by Wiister et al. (2010). 

Free space optical lattices have been formed in many geometries from one dimen- 
sional to three dimensional. In a parallel development also evanescent fields have 
been used to trap atoms and create new lattice geometries. In a particularly interest- 
ing and potentially important example, tapered optical fibers were used to trap neu- 
tral atoms (Bhagwat and Gaeta, 2008; Christensen et al., 2008; Bajcsy et al., 2009). 
These are coupled to the propagating fiber modes via evanescent fields surrounding 
the fiber (Nayak et al., 2007). Such a hybrid combination of atoms and a solid state 
devices allows trapping and, simultaneously, optically interfacing (Hammerer et al., 
2010). Fibers with diameters smaller than the wavelength exhibit strong transverse 
confinement and a pronounced evanescent field, which by using two colors can be 
designed to create an array of deep optical micro-traps. It was recently realized 
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with Cesium atoms, where it was shown that the atoms can be efficiently trapped 
and interrogated with a resonant light field sent through the nanofiber (Vetsch et al., 
2010; Goban et al., 2012). This setup provides an ideal test system to study one di- 
mensional excitons, for which we developed a microscopic theory. The strong col- 
lective coupling between atoms and field modes naturally leads to the formations 
polaritons in this setup (Zoubi and Ritsch, 2010c). As one of the most interesting 
consequences we could show that one can create polaritons with a strong fiber mode 
contributions but with virtually no coupling to the free space radiation modes, this 
should constitute an ideal photon-atom ensemble interface for coherent quantum 
information storage and processing. 

The review is organized as follows. In section 2 we extend the Bose-Hubbard model 
to exploit electronically excited atoms. In section 3 we introduce excitons into a 
system of ultracold atoms in optical lattices. We discuss their formation, dispersion 
and life times for different geometries and dimensionality. Section 4 is opened by 
discussing the quantum phase transition within a cavity followed by introducing 
cavity polaritons in the strong coupling regime, where various possible physical 
realizations are studied. Defects of atoms excited into higher Bloch bands or va- 
cancies in the Mott insulator phase are treated in section 5. The conclusions are 
contained in section 6. 



2 Ultracold Atoms in an Optical Lattice as Artificial Crystals 

2.1 Superfluid to Mott-Insulator Transitions 

Optical lattices created by off resonance laser fields constitute a regular array of 
micro-traps into which a dilute gas of ultracold atoms can be loaded (Greiner et al., 
2002; Bloch et al., 2008). The atoms experience an optical lattice potential with lat- 
tice constant of half the wave length of the lattice laser with different depth for each 
internal atomic state. At ultra low temperatures, T « 0, atoms are localized at the 
lowest Bloch band of the lattice. In such structures, atoms can hop among different 
sites through quantum tunneling and interact locally via scattering. The tunneling 
process is mainly among nearest neighbor sites with the hopping parameter J, and 
the atom-atom interactions are efficient only for atoms at the same lattice site with 
an effective strength U, which is taken to be a repulsive determined by the s-wave 
scattering length. This system is well described by a Bose-Hubbard model (Jaksch 
et al., 1998) represented by the Hamiltonian 
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Fig. 1. Phase diagram: (ji/zJ) vs. (U/zJ). The Superfluid SF and the three Mott insulator 
MI regions, for n = 1,2, 3, are shown. Figure reprinted with permission from Zoubi and 
Ritsch (2009d) Copyright 2009 by American Physical Society. 



where b] and bi are the creation and annihilation operators of an atom at site i, and 
s t is the atom on-site energy, which includes the external trap potential, n is the 
chemical potential accounting for a possible exchange of atoms with an external 
atomic reservoir. The i summation is over the lattice sites, and the bracket indicates 
a possible hopping only between nearest neighbor sites. 



The Bose-Hubbard model predicts a quantum phase transition between the super- 
fluid and the Mott insulator phase by changing the external laser field intensity, that 
is changing the kinetic versus the repulsion energy of the quantum gas. In the limit 
of J » U the atom hopping among the lattice sites dominates the dynamics and 
the atomic states are spread out over the whole lattice and one gets a superfluid. 
In the opposite limit J «; U, the on-site atom-atom interaction dominates and the 
atoms localize deterministically at individual sites, so that we have a Mott insula- 
tor. The corresponding phase diagram can be calculated in the mean-field approach 
(Sheshadri et al., 1993; Sachdev, 1999; van Oosten et al., 2001). In Fig. 1 the phase 
diagram is plotted in the plane (/i/z/) - (U/zJ) which is scaled by zJ, where z is the 
number of nearest neighbors. The Superfluid (SF) phase appears outside the shaded 
three Mott Insulator (MI) phase regions, which are for one, two, and three atoms 
per site, that is n = 1,2,3, where we assumed e, = 0. In the mean field approx- 
imation the systems dimensionality enters only through z, where in a hyper-cubic 
lattice for the one dimension d = 1 we have z = 2, for d = 2 we have z = 4, and for 
rf = 3we have z = 6. 
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Ground state optical lattice 



Fig. 2. Schematic plot of the ground and excited state optical lattice potentials. 
2.2 Mott Insulator for a Two-Component Bose-Hubbard Model 

In the present review we study collective electronic excitations of an atomic Mott 
insulator. Hence we have to introduce an excited atom state, which modifies the 
above phase diagram as we now have two-level atoms treated here as two kinds 
of atoms, namely ground state atoms and atoms excited to a fixed internal state. 
When loaded on an optical lattice due to different polarizability and Stark shift of 
the ground and excited states each group of atoms experiences a different optical 
lattice potential, but with ground and excited potentials showing extrema at the 
same positions, as schematized in Fig. 2. For certain, magic, laser wave lengths, 
light shift yield nevertheless identical potentials. This is important e.g. to create an 
optical lattice clock (Takamoto et al., 2005). Close to T = the ground and the 
excited atoms will again occupy only their lowest Bloch band. 

We consider the two groups of atoms with different internal states, hyperfine or 
electronic, as two kinds of bosons. Therefore, we generalize the single-component 
Bose-Hubbard model into two-component Bose Hubbard model (Jaksch et al., 
1998; Chen and Wu, 2003; Zoubi and Ritsch, 2009d), with the Hamiltonian 



H = -J g J] b]bj - /. £ cjcj + £(ef - n g ) b% + J>« - He) c\ Ci 

(i,J) <Uj) i i 

+ y £ b ' b ' bibi + T Z c " c ' CiCi + u « Z b * c i biCi > (2) 

i i i 

where b\ and b t are the creation and annihilation operators of a ground state atom 
at site i, respectively, with the on-site energy e g t , and chemical potential fi g ; c] and 
Ci are the creation and annihilation operators of an excited state atom at site i, 
respectively, with on-site energy s e r and chemical potential fi e ; where = ef + 
hco a , with oj a being the effective atomic transition frequency, which can include AC 
Stark-shifts relative to the free atomic transition frequency. Here J g and J e are the 
hopping parameters of ground and excited state atoms, respectively. U g and U e are 
the on-site ground and excited state atom-atom interactions, respectively, and U eg 
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Fig. 3. Phase diagram: (jj. g /zJ g ) vs. (U g /zJ g ), and (jj. e /zJ e ) vs. (U e /zJ e )- We used 
U eg lzJg = U eg lzJ e = 15, with e? /zJ g - and ef/z/ e = 100, for ra g = n e = 1. Beside 
the SF and the MI regions, the SM region is associated to the Mott insulator of excited 
atoms and superfluid phase of ground atoms, and the MS region to the Mott insulator of 
ground atoms and superfluid phase of excited atoms. The full-line Lg determines the SF-MI 
transition of the ground atoms, and the dashed-line Le that transition for the excited atoms. 
Figure reprinted with permission from Zoubi and Ritsch (2009d) Copyright 2009 by Amer- 
ican Physical Society. 



is the on-site ground-excited atom interaction related to scattering between ground 
and excited state atoms. This extended model leads to a much richer phase diagram. 



Now we present one typical example of the phase diagram, which is calculated in 
using the mean-field approach. In figures 3 we plot the phase diagram for the plane 
(/Ug/zJg) vs. (Ug/zJg), and (jj. e /zJ e ) vs. (U e /zJ e ), which are scaled by zJ g and zJ e - 
We used for the energies s^/zJ g = and s'j/zJe = 100, and for the scaled atom-atom 
coupling the number U eg lzJ g = U eg /zJ e = 15. Here the excited and ground phase 
transition lines split, the transition line (Le) is for excited state atoms, and (Lg) 
for ground state atoms. Beside the superfluid phase, we get three other regions. 
The Mott insulator region with one ground state and one excited state atom per 
site, that is n g = n e = 1. The SM region, where the excited state atoms are in 
the Mott insulator phase with one atom per site, and the ground state atoms are 
in the superfluid phase. The MS region, where the ground state atoms are in the 
Mott insulator phase with one atom per site, and the excited state atoms are in the 
superfluid phase. We like to mention here that more complex quantum phases are 
suppressed under the present mean field theory, e.g., that of the anti-ferromagnetic 
type in which the Mott insular includes an excited atom at one site and a ground 
one in the nearest neighbor sites (Altman et al., 2003). 
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3 Excitons in Optical Lattices 

As shown above even for two state atoms one can prepare a perfectly filled Mott 
phase for the atoms, which we will use as the basis of our further considerations 
involving resonant excitations. We focus on a single atom per site, but eventually 
allow for two atoms per site as in Subsection 3.5. The Mott insulator constitutes 
an artificial crystal with a lattice constant much larger than the average dimen- 
sion of the atom. Therefore, each atom retains its identity except from the internal 
shifts and there are no overlaps among electronic wave functions of atoms setting 
at different sites as in conventional solids. We conclude that our artificial lattice is 
similar to molecular or noble atom crystals but allow for the wide controlablity of 
the system parameters, as the artificial lattice is formed by variable external laser 
field. Solid molecular crystals, however, are fixed by the chemical bonds, e.g. van 
der Waals forces, and can be only slightly controlled in applying external fields. 
Frenkel excitons play a central role in optical and electrical properties of molecular 
crystals (Davydov, 1971; Agranovich, 2009; Zoubi and LaRocca, 2005a). Here we 
ask the question related to the possibility of the formation of Frenkel like-excitons 
for optical lattices. First we discuss the possibility of energy transfer among atoms 
at different sites induced by electrostatic interactions or resonant dipole-dipole in- 
teractions, which is the main mechanism behind the formation of excitons. 



3.1 Resonance Dipole-Dipole Interactions 

Lets take the simplest case of two atoms setting on a lattice at nearest neighbor sites 
such that they are separated by the lattice constant a. The atoms are considered to 
be two-level systems with transition energy E A = hco A . If one atom is excited and 
the second is in the ground state, can this excitation be transferred among the two 
atoms, in which the excited atom decays to the ground state and the ground state one 
becomes excited. Here we consider only energy transfer without charge transfer? 
This process is possible through the resonance dipole-dipole interaction (Craig and 
Thirunamachandran, 1984) with the energy transfer parameter 



where q A is the atomic transition wave number, as E A = hq A , and T A is the single ex- 
cited atom damping rate (Loudon, 2000), which is given by Y A = a) 3 A fi 2 (3neohc 3 )~ l , 
with n the transition dipole. Here 6 is the angle between the transition dipole and 
vector distance between the two atoms. 




(3) 
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Fig. 4. The scaled interaction J (qAa) /Ia vs. qAa, for (a) 9-0° and (b) 9 = 90°. The 
full line is for equation (3), and the dashed line for equation (4). Figure reprinted with 
permission from Zoubi (2012a) Copyright 2012 by European Physical Society. 

The energy transfer is efficient in the limit of \J (^a) \/Y A > 1, where the electronic 
excitation can transfer among the atoms before it spontaneously decays, and this 
limit is achievable up to q A a ~ 1. For a typical optical lattice we have, e.g., Ea = 
1 eV, with q A * 4 x 1(T 4 A -1 . For a = 1000 A we get q A a ~ 0.5. For 9 = 0° we 
obtain J(0.5)/T A * -13.4, and for 6 = 90° we obtain 7(0.5)/r A * 5.4, in which 
the energy transfer is possible. Note that in the limit q A a < 1, with efficient energy 
transfer, the electrostatic interaction limit is applicable (Zoubi, 2012a). In this limit 
and in neglecting the radiative corrections, we get 



Using the previous numbers, 6 = 0° yields J(0.5)/T A ~ -12, and 6 = 90° yields 
J(0.5)/Y A ~ 6, which are slightly different from the above exact results. For smaller 
qAa we get much better agreement. Figure 4(a) reports J (qAd) IT a as a function 
of q A a for Eqs. (3) and (4) at the polarization direction = 0°, and in (b) for 
6 = 90°. The results justify the use of electrostatic dipole-dipole interactions for 
optical lattices when q A a < 1. 

Energy transfer of electronic excitations among close lattice sites induced by dipole- 
dipole interaction tends to delocalize the excitation coherently in the lattice. Ex- 
ploiting the lattice symmetry we can construct collective excitations in such sys- 
tems, which are called excitons. In momentum space such de-localized electronic 
excitations can be represented by propagating excitation waves. In the following 
we explicitly calculate these excitons for various optical lattices with different ge- 
ometries and dimensionality. 

We concentrate now in the Mott insulator phase with one atom per site (Zoubi and 
Ritsch, 2007). The atoms as before are assumed to be two-level systems, and the 
electronic excitation can transfer among atoms at different sites, e.g. among sites 
n and m, with the transfer parameter J nm . As discussed above, such a process can 
be efficient in the limit of \ J nm \ > Y A , where an electronic excitation could transfer 




(4) 
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before the decay of the excited atom. The excitation Hamiltonian in general can be 
written in the second quantized form as 



H = Y j E a B% + 



^ %J nm B f n B, 



(5) 



;? 



nm 



where B n and B n are the creation and annihilation operators of an electronic exci- 
tation at site n. The first term is for the on-site excitation and the second for the 
energy transfer. We neglect terms that do not conserve the number of excitations. 
The operators are of the spin-half type, and they behave on-site as fermions, and 
as bosons among different sites. But for a single excitation, or for very small num- 
ber of excitations relative to the number of lattice sites, saturation effects can be 
neglected, and then the possibility of having two or more excitations to appear at 
the same lattice atom is negligible (Agranovich, 2009; Zoubi and LaRocca, 2005b). 
Therefore, from this point on we assume the excitations to behave as bosons with 
the commutation relation [B n , B} n ] = 6 nm . 

3.2 One -Dimensional Atomic Chains 

We treat first a long one dimensional optical lattice (Zoubi and Ritsch, 2010b), as 
in Fig. 5. The lattice length is enough long to neglect edge effects, which allows 
the use of the periodic boundary condition. The previous Hamiltonian can be easily 
diagonalized in applying the Fourier transform 



where N is the number of lattice sites, and z„ = an is the position of an atom at site 
n. Here k is a wave number that in using periodic boundary conditions takes the 
values k = with (p = 0, +1, • • • , +N/2), and L = Na is the lattice length. This 
solution holds also for a ring of an atomic chain. In applying the transformation the 
diagonal Hamiltonian reads 



where the exciton dispersion is given by E k = E A + J k . As the energy transfer is a 
function of the interatomic distance, R = a(n - m), the exciton dynamical matrix 
is given by E k = }Zr Jk(R)e' kR - Note that in place of a discrete energy level of the 
excited atom we get an energy band for the excitons. 

The single exciton state of wave number k is obtained from the vacuum state, |0) = 
\{g„}) in which all atoms are in the ground state, in applying the exciton creation 
operator, namely \l k ) = B'[\{g n )), where \l k ) = -j= \e n )e~ ikan . Here \e n ) stands for 
a state where only atom n is excited and all the others are in the ground state. The 




(6) 




(7) 



k 
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Fig. 5. Schematic plot of a one dimensional atomic chain of two-level atoms of transition 
energy E A , lattice constant a, and energy transfer parameter J. 




Fig. 6. A one-dimensional infinite lattice of lattice constant a, and with one atom per site. 
The transition dipole // is plotted which makes an angle G with the lattice axis. For a finite 
optical lattice of N sites with one atom per site, the lattice length is L - a(N +1), 

single exciton is a coherent state built of a superposition of electronic excitations 
at all sites, where the probabilities to find an excitation at each site are equal. This 
state is a collective electronic excitation and is considered as a wave that propagates 
to the left or the right in the lattice with a fixed wave number k. A superposition of 
such states can be used to build a more physical state of a wave packet of excitons. 

We assume now energy transfer only among nearest neighbor sites. In the q A a < 1 
limit, the transfer parameter 7 is given by Eq. (4) and shows a strong dependence 
on the 6 angle between the directions of the transition dipole and of the lattice. 
The transfer parameter changes from being negative (attractive) at 8 = to positive 
(repulsive) at 9 = n/2. The dispersion relation now reads 

E k = E A + 2Jcos(ka), (8) 

with a bandwidth of 47. For small wave numbers, that is in the limit of ka «; 
1, or for long wavelength excitons, we get the dispersion E k « E A + 27 + j^-, 

where we defined the exciton effective mass m ex = j^. The exciton in this limit 
is a quasi-particle that propagates in the lattice with a given wave number and an 
effective mass. At k « we get E k « E A + 27, with a shift produced by the dipole- 
dipole interaction. In this limit of long wavelength excitons, beyond the nearest 
neighbor interactions the long range dipole-dipole interaction yields the correction 
E k « E A + 2^(3)7, where f(3) * 1.2 (Zoubi and Ritsch, 2010c). Such a shift has 
a large influence on optical lattice clock states (Takamoto et al., 2005), but can be 
controlled through the direction of the transition dipole (Chang et al., 2004), and 
even can be brought to zero at the angle of 6 ~ 54.7. 

Next we treat a finite optical lattice chain of iV two-level atoms (Zoubi and Ritsch, 
2009c, 2012b), as schematized in Fig. 6. The lattice length is taken to be smaller 
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Fig. 7. A planar square optical lattice with lattice constant a. The directions of the in-plane 
wave vector k and the transition dipole // are indicated. 

than the atomic transition wavelength. In considering energy transfer only between 
nearest neighbor atoms, the Hamiltonian can be diagonalized by introducing stand- 
ing wave excitons 



with (k = 1, • • • , N). As a mathematical trick to find appropriate solutions, we add 
two additional empty sites at the two lattice edges, at n = and n = N + 1, where 
the collective electronic excitation wave function vanishes. The exciton dispersion 
then reads 



The single electronic transition energy splits into N discrete collective electronic 
excitation energies. The modes are divided into two groups, symmetric and anti- 
symmetric modes. In the next section we will show that the antisymmetric modes 
are dark, and the symmetric ones are bright. Furthermore the first bright mode is 
node less and found to be super radiant and dominates the optical properties. 



3.3 Two -Dimensional Planar Optical Lattices 

Let us consider the case of a two dimensional planar optical lattice (Zoubi and 
Ritsch, 2007). For simplicity we treat the case of a square symmetry with lattice 
constant a, and one two-level atom per site, as in Fig. (7). The position of a lattice 
site is defined by the in-plane vector n = (n x ,n y ), where n t = a(0, ±1, • • • , +NJ2), 
with N = N x x N y the number of lattice sites. The electronic excitation Hamiltonian 
is given by 




(9) 




(10) 




(11) 



n 



nm 



which can be diagonalized by applying the transformation 




(12) 
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to get the diagonal exciton Hamiltonian 

H« = £)£«(k)4fl k , (13) 

k 

with the energy dispersion 

E ex (k) = E A + Yj We ik R , (14) 

R 

where we defined the distance between two sites by R = n - m and introduced as 
previously the 7 transfer parameter. 

The exciton is now represented by a wave that propagates in the lattice with in-plane 
wave vector k, where we have k { = jf - (0, ±1, • • • , ±Ni/2). The nearest neighbor 
interaction yields the dispersion with E ex (k) = E A + 2J[cos(k x a) + cos(k y a)], which 
for small wave numbers that is ka <sc 1 with k = |k|, reads E k « E A + 47 + 

with the exciton effective mass m ex = At k « we get E^ ~ E A + 47 with 
the dipole-dipole interaction shift. In the limit of long wavelength excitons, beyond 
the nearest neighbor interactions the long range dipole-dipole interaction yields the 
correction Ek ~ E A + IF J, where in using the Ewald's summation method (Zoubi 
and Ritsch, 201 lb) we get F * 9/2. 

3.4 Radiative Decay of Excitons 

The exciton damping rate through the spontaneous emission into free space can be 
calculated in using Fermi's golden rule (Zoubi and Ritsch, 2010b, 2011a, 2012b; 
Zoubi, 2012a; Ostermann et al., 2012) and shows significant deviation from a single 
excited atom (Loudon, 2000). For a finite atomic chain the standing wave excitons 
are divided into two groups of dark and bright modes, which exhibit metastable or 
superradiant propagating excitons. Let us mention here an important point concern- 
ing the atomic recoil effect via the emission of a photon. As the exciton is a collec- 
tive electronic excitation delocalized in the lattice, where the exciton momentum 
is shared among all lattice atoms, the emitted photon inherits the momentum and 
energy of the exciton, which requires the photon wave vector component parallel 
to the lattice to be equal to the exciton wave vector. In this process the whole lattice 
reacts and absorbs the momentum collectively and not only a single atom. Hence, 
the recoil heating is significantly suppressed, and we can conclude that excitonic 
interactions will not heat the lattice significantly and destroy the Mott insulator. 

We start with a long optical lattice with one atom per site (Zoubi and Ritsch, 2010b) 
as in Fig 6. We assume that a single exciton is excited in the lattice with a fixed wave 
number k, and a given polarization 9 for the transition dipole direction relative to 
the lattice axis. The damping rate of the exciton into the radiation field of the free 
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Fig. 8. The damping rate vs. ka, for (a) 9 - and (b) 9 = 90. The damping rate is zero 
beyond the critical wave number k c a ~ 0.507. The dashed-line is for a single excited atom 
damping rate. Figure reprinted with permission from Zoubi and Ritsch (2010b) Copyright 
2010 by European Physical Society. 

space is given by 

H 2 E 2 k { , {tick) 2 , , 

where Ek is the previous one dimensional exciton dispersion, e.g. as for the nearest 
neighbor interactions in Eq. (8). For some wave numbers k and certain angles 9 
the damping rate can be much larger than the single excited atom one, as we get 
a superradiant enhancement analogous to the Dicke model. However, for other k-s 
and 9-s the damping rate can be much smaller than the single excited atom one, or 
even completely vanish. In general the damping expression shows the existence of 
a critical wave vector k c that is given by E ex (k c ) = tick c , beyond which the damping 
rate is zero and the excitations can no longer decay radiatively. For 6 = 0° the 
damping rate at the critical wave vector k c is zero, and for 6 = 90" the damping 
rate at k c is finite and equal to T(k c ) = /j. 2 E 2 x (k c )(2eQati 3 c 2 y l . We conclude that 
Eq. (15) holds up to the critical wave number, k c , where beyond k c no free space 
photon exists which conserves both energy and momentum component parallel to 
the lattice. 

We exhibit this behavior in two plots of the damping rate using the following pa- 
rameters: transition energy E A = 1 eV, lattice constant a = 1000 A, and transition 
dipole [i = 1 eA. In Fig. 8(a) we plot as a function of ka for 9 = 0° the damping 
rate compared to the damping rate of a single atom. Note that for small wave num- 
bers, ka ~ 0, the exciton damping rate is much larger than the free atom damping 
rate, and these states are superradiant states. With increasing wave numbers the 
damping rate drops and reaches the single atom rate, and beyond a critical wave 
number k c the damping rate becomes zero and no damping is obtained beyond k c . 
The critical wave number is obtained here at k c a « 0.507. In Fig. 8(b) the same 
plot is for 9 = 90°. Here the damping rates are larger than the single atom one for 
small wave numbers, ka ~ 0, but now the damping rate increases with increasing 
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Fig. 9. The r^/T^ collective electronic excitation damping rate rapported to that of a sin- 
gle excited atom k, for N = 10. Figure reprinted with permission from Zoubi and Ritsch 
(2012b) Copyright 2012 by European Physical Society. 



the wave number. Namely, the large wave number states became more superradiant 
states. The maximum damping rate is obtained at the critical wave number k c , and 
the damping rate is plotted up to that value, as beyond this point the excitons are 
metastable. Note that excitons can appear up to the boundary of the Brillouin zone 
at k B a = n, and hence excitons in the region k c < k < k B are metastable. These 
excitons are not dark, only they are decoupled from the free space photons, but 
they can be coupled to confined photons, as presented in Subsection 4.5 through 
the coupling to tapered nanofiber photons (Zoubi and Ritsch, 2010c). 

The results are drastically changed for a finite one dimensional atomic chain of 
only a few atoms (Zoubi and Ritsch, 2009c, 2012b; Zoubi, 2012a; Ostermann 
et al., 2012) which is short relative to the atomic transition wave length. In this 
case we obtained that the exciton states are divided into two groups of symmetric 
and antisymmetric standing waves. For each collective mode we can define a col- 
lective transition dipole, which is assumed to be localized at the center of the lattice, 
for the limit of short chains. The antisymmetric excitons with even wave numbers 
(k = 2,4, • • •) are dark with zero damping rates, that is Y k = 0. The symmetric 
excitons with odd wave numbers (k = 1, 3, • • • ) are bright with the damping rates 

lx 2 E\ (2a\ Jnka\ 

where the lattice length is L = a(N +1), and the exciton dispersion is given in 
Eq. (10). The first mode (k = 1), the one without nodes, is superradiant, with a 
damping rate nine times larger than the second bright mode (k = 3). In Fig. 9 the 
collective electronic excitation damping rate rapported to that of a single excited 
atom, T k /T A is plotted versus k, using the previous parameters, 6 = and N = 10. 
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Interesting results can be also obtained for planar optical lattices (Zoubi and Ritsch, 
2011a). We deal here, as before, with two dimensional large optical lattice with 
one atom per site and of square symmetry, as plotted schematically in Fig. 7. We 
calculate the damping rate into free space of an exciton with in-plane wave vector 
k. Now the dipole makes an angle 6 with the normal to the lattice plane, and an 
additional angle <p is needed that is between the in-plane component of the dipole 
and k. The result is 



r.(k) = -4 - E JM W^-cos^) 

2e ° a n c ^(k) - Efck) I V 

+ cos2 d 15T7n - 2 sin cos cos J^- J^Jk) - E 2 Q (k)\ , (17) 
E 2 ex (k) E 2 ex (k) y J 

where the exciton dispersion is defined in Eq. (14), and E Q (k) = Tick. The most 
important result here is that for k > k c we have Y k = 0, where tick c = E ex (k. c ), 
as beyond the singularity point the damping rate becomes imaginary. The results 
should be compared with the radiative damping rate of a single excited atom r^. 

We present numerical results using the typical previous numbers for optical lattices, 
except for simplicity we neglect the k dependence of the exciton energy, negligible 
compared to the E (k) one. Fig. 10(a) shows T ex /T A vs. E (k) for 9 = 0, where 
the transition dipole has only component normal to the lattice plane. Here T ex (k) 
is tp independent. The damping rate starts from zero, where the long wave length 
excitons are metastable, and increases with k to become superradiant with damping 
rate larger than T A . Close to fick c = 1 eV the damping rate diverges, and beyond 
that value it jumps back to zero, all excitons between k c and the Brillouin zone 
boundary at n/a being metastable. Fig. 10(b) shows T ex /T A vs. Eo(k) for 6 = n/4, 
where the transition dipole has equal in-plane and normal components. We take 
the in-plane wave vector k to be parallel to the in-plane dipole, that is (f> = 0. 
Very interesting behavior appears here for the damping rate. For small in-plane 
wave numbers the excitons are superradiant with damping rate larger than Y A , and 
decreases in increasing k till it becomes zero at E (k) = E ex (k)/ V2 where excitons 
around this point are metastable. For larger k the damping rate start to increase and 
diverges at E (k c ), and then jump back to zero beyond this point. 

The dark and metastable states have an important potential use for the quantum 
information processing and the memory devices (Bouwmeester et al., 2000) but 
their excitation is an experimental challenge. The dark state can be addressed by 
applying external fields or by adding a structural asymmetry. The metastable ones 
can be coupled to confined photons in using different optical resonators as shown in 
Subsection 4.5. The bright states are coupled directly to photons through the dipole 
coupling, and due to their strong dipole moment they can easily be tuned to achieve 
the strong coupling regime, making the atomic chain an optimal candidate as an 
active material for cavity QED. 
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Fig. 10. The scaled damping rate T ex /TA vs. Eo(k) = Tick, for (a) 6 = 0, which is <f> in- 
dependent, and for (b) 6 = n/4 with = 0. Here ftcn/a ~ 6.2 eV. Figure reprinted with 
permission from Zoubi and Ritsch (2011a) Copyright 2011 by American Physical Society. 

3.5 Excitons in Optical Lattices with Two-Atoms per site 



Up to this point we have treated only the case of an optical lattice in the Mott in- 
sulator phase with one atom per site, but the case of two atoms per site is of large 
interest and can be easily realized in the present optical lattice experiments (Spiel- 
man et al., 2007). Furthermore, this case shows novel phenomena beyond that of 
one atom per site, e.g. the appearance of dark and localized excitons (Zoubi and 
Ritsch, 2008b). As before we consider two-level atoms, in which the ground and 
excited optical lattice potentials have minima at the same positions, and the atoms 
are located in the lowest Bloch band. Moreover, we neglect overlaps between elec- 
tronic wave functions for on-site and off-site atoms. Namely, we neglect the on-site 
formation of molecules from the two atoms, and we consider only transfer of elec- 
tronic excitations. For planar optical lattices, the electronic excitation Hamiltonian 
is as in Eq. (11), but with an additional index for the two on-site atoms 

H ex = j]e a b:'b: + « (i8) 

where a and B run over the two on-site atoms, which are denoted by (1) and (2). 
For simplicity, we assume energy transfer only among on-site atoms with average 
coupling parameter Jj^ = J^ a = J , and among nearest- site atoms with coupling 
parameter /^ m) = 7^ m) = y^ m) = 7 < 2 n 1 m) = J, and we consider a lattice of square 
symmetry. The above Hamiltonian can be diagonalized in two steps, first the diago- 
nalization related to the two on-site atoms, and second related to the site positions. 
The first transformation is taken in introducing the on-site symmetric and antisym- 
metric states 

B s n = S— i K = S— i, (19) 
V2 V2 

to get 

H„ = J] {E* B«B s n + E a B a :B a D ) + J]2J B«B* m , (20) 

n (nm) 
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Fig. 11. (a) The two-level on-site two atoms with a transfer parameter Jq, and a transition 
energy E A . (b) The two diagonal states with the splitting energy of 2Jq. 

with the energies E s = E A + J and E a = E A - Jq. The on-site energy states are 
represented in Fig. 11. The second transformation is given by B v a = -4= £k e lk n B^, 
which yields the diagonal Hamiltonian 

H ex = Yj ^v(k) flfflj + Yj Ea B » B »- (21) 

k n 

The on-site symmetric states are delocalized in the lattice and form excitons with 
energy dispersion E s (k) = E A + Jq + 4j[cos(k x a) + cos(& y a)j. The antisymmetric 
states are on-site localized and form no excitons in the lattice; they are dispersion- 
less with energy E a = E A - J . The on-site symmetric state is found to be super- 
radiant with damping rate of two time a single excited atom, that is Y s = 2T A , and 
the antisymmetric state is dark with Y a = 0. These results can be obtained from 
Eq. (16) of Subsection 3.4 for a finite chain with N = 2, as presented also in Ficek 
and Tanas (2002). 

Dark localized states can in principle also be excited through nonlinear optical 
processes or the decay of multiply excited atomic states (Ostermann et al., 2012). 
Moreover, the addition of asymmetry to the system through the application of ex- 
ternal fields or the deformation of the lattices sites, e.g. in using different laser fre- 
quencies and intensities for different directions, can make the antisymmetric states 
to be slightly bright and to delocalize in the lattice. This asymmetry can be applied 
in a controllable way for the writing and reading of the antisymmetric states. 

Another possibility is the use of an optical super-lattice with two on-site potential 
minima. A single atom is taken to be localized at each minimum, and the two on-site 
atoms are now separated by a given distance. This leads to the exciton Hamiltonian 
Hex = Zkv£y(k) ^k' where (v = a, s) stands for the two exciton branches, with 
the symmetric state dispersion E s (k) = E A + J + AJ S [cos(k x a) + cos(^a)j, where 
J S = J' + J", here 7 < 1 n 1 m> = 7f n 2 m> = J' and j£ m) = J% m) = J". The antisymmetric 

dispersion E a (k) = E A - J Q + 4J a [cos(k x a) + cos(^a)j, where J a = J' - J". The 
antisymmetric state can transfer among nearest neighbor sites with the parameter 
J a , and to be delocalized and to form a second exciton branch, with a narrow energy 
band, and is slightly bright with a finite damping rate. 

The formation of an on-site molecule out of the two on-site atoms significantly 
changes the optical properties, as the diatomic molecules can be oriented on the 
lattice due to the laser polarization or the external trapping fields. This makes the 
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Fig. 12. Two parallel one-dimensional optical lattices which are separated by a distance d, 
with one atom per site and lattice constant a. The transition dipole makes an angle 9 with 
the lattice direction. 

system anisotropic similar to uni-axial crystals (Zoubi and Ritsch, 2009b). As pre- 
sented with more details in Subsection 4.3, the localization of the anisotropic lattice 
between planar cavity mirrors results in cavity photon polarization mixing, which 
can be observed through linear spectra. 



3.6 Coherent Transfer ofExcitons among Optical Lattices 



In this Subsection we discuss the exciton coherent transfer between two different 
optical lattices (Zoubi and Ritsch, 2010c, 2011b, 2012b). We start with two atoms 
that are separated by a vector distance R and have transition dipoles n a and «g. 
The transition dipoles are equal for identical atoms, but they can have in general 
different orientations in space. The dipole-dipole interaction is given by 

^ (R) = 4^;l^ w r (22) 

which behaves as inverse cube of the inter-atomic distance. This electrostatic in- 
teraction holds for inter-atomic distances smaller than the atomic transition wave 
length A A , where Ea = Hc/Aa- This interaction can be used also for two finite atomic 
chains, with lengths and inter-chain distances smaller than A A , where the transition 
dipoles are taken to be of the collective electronic excitations. The energy transfer 
can be blocked among chains of different lengths as the collective state energies 
can be off resonance (Zoubi and Ritsch, 2012b). 

Because things can be completely different for extended systems, first we present 
the coherent exciton transfer among two long parallel one dimensional optical lat- 
tices (Zoubi and Ritsch, 2010c), as schematized in Fig. 12. The two lattices are 
separated by a distance d, with one atom per site and lattice constant a. We con- 
sider an exciton with a given wave number k that propagates in one optical lattice 
with a fixed polarization 6, and we derive the coupling parameter for the exciton 
transfer into the other lattice. We calculate this parameter in considering the long 
range dipole-dipole interaction among each pair of atoms in applying the Ewald's 
summation method of Born and Huang (1954). In the limit of small wave numbers, 
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Fig. 13. Two parallel planar optical lattices that are separated by a distance b, with one 
atom per site and lattice constant a. The transition dipole includes in-plane and normal 
components. 

that is ka «; 1, and for large inter-lattices distance, that is kd » 1 so that d » a, 
we obtain 

J(k, 6) = e~ k V [2 cos 2 6 - 1} . (23) 

The formation of excitons affect strongly the transfer parameter, as in place of the 
inverse cube dependence on the inter-atomic distance for the transfer parameter we 
get exponentially decay parameter behavior. 

Now we consider the case of two planar parallel optical lattices (Zoubi and Ritsch, 
2011b), which are separated by a distance b, as schematized in Fig. 13. For an 
exciton with in-plane wave vector k that propagates in one lattice, we calculate 
the transfer parameter into the other lattice in applying once again the Ewald's 
summation method. The coupling parameter is 

J '^ = ^^{{^^f -Ml}kae- kb , (24) 

where we defined // = (ji\\,fx z ) with the in-plane dipole fi^ = (pi x ,fi y ), and the unit 
vector k = k/fc. This result shows a strong anisotropic effect. 

The most important feature here is the exponential decay behavior as a function of 
the inter-lattices distance of the exciton coupling among the two optical lattices. 
Even though the resonance dipole-dipole interaction has inverse cubic dependence 
of the distance between two localized dipoles, the appearance of excitons induce 
exponentially decrease dependence for the exciton collective dipole among the two 
lattices. For several parallel lattices the result allows the use of the assumption of 
energy transfer among only nearest neighbor planes. 

Such coherent coupling of the exciton transfer among the two lattices can be re- 
sponsible for the formation of hybrid excitons in the whole systems, in which an 
exciton can coherently oscillate between the two lattices. The excitation Hamil- 
tonian for two parallel planar optical lattices is exactly as in Eq. (18), but here a 
and 8 stand for the two lattices. The same diagonalization steps lead to the exci- 
ton Hamiltonian H ex = £kv E v (k) B^B^, with two hybrid exciton branches. The 
symmetric excitons have the dispersion E s (k) = E A + 7(k) + 7'(k), and the anti- 
symmetric dispersion is E a (k) = E A + 7(k) - 7'(k). Here 7(k) is the intra-lattice 
coupling, and /'(k) is for the inter-lattices, which is given in Eq. (24). Note that 
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the two hybrid excitons have an energy split of 2J'(k). Also here the antisymmetric 
modes are dark with zero damping rate, and the superradiant ones are bright with 
two times the single exciton damping rate. 



4 Cavity QED with Excitons: Polaritons 

In the recent years few experiments implemented a dilute boson gas of ultracold 
atoms in an optical cavity (Brennecke et al., 2007; Colombe et al., 2007; Slama 
et al., 2007), and they reached the strong coupling regime because the collective 
coupling enhancement is enhanced by the square root of the atom number. In the 
experiment of Ritter et al. (2009) the optical lattice was formed by off resonance 
cavity photons, with the resulting force effects investigated in Maschler and Ritsch 
(2005); Mekhov et al. (2007); Larson et al. (2008). Optical lattices within an op- 
tical cavity operating in the resonant regime are not yet implemented. Zoubi and 
Ritsch (2007, 2008b, 2009c) have introduced optical lattices as a new type of active 
materials for cavity QED. Beside the importance of such a system for fundamental 
physics and its wide range of applications, the appearance of superradiant collec- 
tive states in optical lattices implies the confinement of the electromagnetic field 
within a cavity. Moreover, as large one and two dimensional lattices give rise for 
metastable states, which decouple to the free space photons, they appeal to confined 
photons for efficient coupling. 



4.1 Quantum Phases for an Optical Lattice within a Cavity 

We begin by answering the crucial question concerning the influence of the cavity 
photons on the previously calculated phase diagrams (Zoubi and Ritsch, 2009d). 
Namely, how the different Mott insulator phases are affected due to the coupling 
to cavity photons. The optical lattice is taken to be localized between cavity mir- 
rors, as in Fig. 14. Two electromagnetic fields are required, the classical laser field 
which is off-resonance and used to form the optical lattice, and the quantized cav- 
ity field. A single cavity mode is assumed to be close to resonance with a single 
atomic transition. The previous two-component Bose-Hubbard model is generalize 
to include the coupling to cavity photons. This coupling changes the behavior of the 
bosons through the emission and absorption of cavity photons, and that is expected 
to modify the phase diagram. 

The two-component Bose-Hubbard model within a cavity in the resonant regime is 
represented by the Hamiltonian 
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Fig. 14. An optical lattice localized along the waist of the Gaussian mode of a cavity with 
spherical mirrors. The lattice length L is smaller than the mode waist Wq. 



where here n g = p g - s g and // e = p. e - s e , with p, g and p, e are the ground and 
excited state atoms chemical potentials, respectively, e? and e\ are taken to be site 
independent, hence we dropped their index i. We use also /i c = p. c - e c , where s c is 
the cavity mode energy and Ji c is the photon chemical potential. Furthermore, the 
photon-excitation coupling / is assumed to be site independent, which is the case 
for a homogeneous optical lattice parallel to the cavity mirrors. We also defined the 
mean atom and photon number operators, by n\ = b]bi, = c]c h and rf = cfia. 
Here a 1 " and a are the creation and annihilation operators of a cavity mode, and 
close to resonance one has e c ~ e e - e g . 

The modified phase diagram is calculated in applying the consistent mean-field 
theory suggested in Sheshadri et al. (1993), as the conventional Bogoliubov ap- 
proximation does not predict the correct quantum phase transition (van Oosten 
et al., 2001). In this mean-field theory the atoms are subject to the mean-field of 
the neighboring sites and of the other kind of bosons. Hence the Hamiltonian can 
be separated into on-site terms and the rest, which can be written as a sum of a 
free part and an interacting part, and where we drop the site index in the follow- 
ing. Then the interaction part can be treated in perturbation theory. In this method 
atom-atom interactions are treated exactly while the hopping and the coupling to 
photons are treated approximately. More complex quantum phases are discussed in 
Bhaseen et al. (2009) beyond the mean-field theory. 



h =- j sYj b ) bi ~ Je Yj c h ~ m§ n < ~ Me ^ 
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We present an example of phase diagram by plotting the scaled chemical potentials 
as a function of the atom-atom interactions. The system in total has a large number 
of open parameters, so that we cannot represent the whole phase space. For our 
arguments here we choose the parameters in such a way to allow the existence 
the Mott insulator phase and to emphasize the modifications brought about by the 
coupling to the cavity photons. We concentrate on the case with two atoms per site, 
that is in the case n e = 1, n g = 1. We assume one photon in average per site, that is 
h c = 1, where h c = n c /N, although this assumption is not so crucial. For n c = 1 we 
have a large number of photons in the cavity, so that photon number fluctuations 
are suppressed. 

As we concentrate in the following on the influence of the coupling to the cavity 
photons, we assume J g = J e = J, and all the parameters will be scaled by zJ ■ 
Furthermore, we assume zero ground state energy, e g = 0, and the cavity mode in 
resonance with the electronic excitation, that is s e = s c . Moreover, by assuming 
a resonant pump-cavity mode the chemical potential for cavity photons becomes 
zero, that is fi c = 0. In the following phase diagram plots we use s c = 100 and for 
the excitation-photon scaled coupling we used F = 25, where F = \f\ 2 h c (z 2 J g Je)~ l - 

Fig. 15(a) reports the phase diagram for the planes fx g vs. U g , and jj, e vs. U e . We 
used U eg = 15 for the scaled atom-atom coupling. Similar to the case without a 
cavity, the plot shows the superfluid and the Mott insulator phases for one ground 
or excited state atom per site, that isn g = n e = 1, with the Le and Lg transition lines . 
Beside the SF and MI Mott regions, we have as before the SM and MS regions. The 
results of Fig. 15(a) should be compared to those of Fig. 3 for the phase diagram 
of a two-component Bose-Hubbard model without coupling to cavity photons. We 
deduce that the effect of the coupling to cavity photons is to shift the Mott insulator 
phase to larger atom-atom interactions, by one order of magnitude in the present 
case. As we used ju c = 0, the ground state transition line is shifted relative to the 
excited state one. 

The coupling to cavity photons tend to delocalize the excited atoms through the 
process in which an excited atom jumps to the ground state by emitting a cavity 
photon, in turn absorbed by another atom in the lattice located farther apart. In 
this situation the atom-atom interaction competes against the atom hopping and the 
coupling to cavity photons, and to recover the Mott insulator phase within a cavity 
one need to apply a deeper optical lattice. 

Fig. 15(b) reports the scaled chemical potentials \x g and fi e as a function of F, for 
U eg = 15, and U g = U e = 250. It is clear that the Mott insulator phase appears only 
for a limited range of coupling parameters, and for a limited range of mean cavity 
photon number per site, which can range from few cavity photons up to about one 
photon per lattice site. Beyond some critical values of F the atom-atom interac- 
tion is not large enough to form the Mott insulator phase and only the superfluid 
phase persists. More Mott insulator phases with other number of atoms per site are 
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Fig. 15. (a) Scaled phase diagram: /j. g vs. U g , and fi e vs. U e , for n g = n e = 1. We have 
£ c = 100, F - 25, U eg = 15, and n c = 1. The dashed line is for the excited state atoms, 
and the full line is for the ground state atoms, (b) Scaled phase diagram: (ju g ,fi e ) vs. F, 
for U g - U e = 250. Figure reprinted with permission from Zoubi and Ritsch (2009d) 
Copyright 2009 by American Physical Society. 

presented in Zoubi and Ritsch (2009d). 



In all the above discussion we assumed long life times for both the excited atoms 
and the cavity photons. This assumption can be a good approximation, e.g., for 
Rydberg atoms in a high-Q cavity. The overlaps between the Mott insulator phases 
of different number of atoms per site make it possible to choose the quantum phase 
such that following the decay of an excited atom the systems may remain in one 
of the Mott insulator phases. The lifetime of the cavity photons has less influence 
on the Mott insulator phase, as the leak of a photon in the Mott insulator keeps the 
system inside the same phase zone. 

The calculated phase diagrams provide us with a new tool for controlling the sys- 
tem quantum phases and their transitions. In place of controlling the external laser 
fields, which form the optical lattice, it is possible to control the cavity photons. 
For example we can chose the different quantum phases by fixing the average num- 
ber of photons within the cavity, which can be easily achieved via controlling the 
external pump. 



4.2 Cavity Polaritons of Planar Two -Dimensional Optical Lattices 



Here we introduce optical lattice ultracold atoms into an optical resonator as an 
active material in a cavity QED (Zoubi and Ritsch, 2007). To represent the main 
results of this set-up we consider a planar optical lattice that enclosed in the middle 
between planar cavity mirrors, as represented in Fig. 16(a). We start by assuming 
an ideal cavity and later we include the damping of cavity photons into the outside 
radiation field. The cavity photons are free to propagate in the cavity plane with a 
given in-plane wave vector k, and they have discrete modes in the perpendicular 
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(b) 




Fig. 16. (a) Cold atom optical lattice located within a cavity of two parallel planar mir- 
rors. The two-dimensional optical lattice created by external lasers is located in the middle 
between, and within a plane parallel to, the two cavity mirrors. The incident, transmitted, 
and reflected fields, with different polarisations are indicated, (b) Schematic plot of an op- 
tical lattice plane with oriented induced transition dipoles. The unitary vectors defining the 
photon polarization, e^, and its orthogonal direction, %, are reported. 

direction with wave numbers k z = mn/L, where L is the distance between the mir- 
rors and m is an integer (m = 1, 2, • • • ). In the perpendicular direction we consider 
only a single mode, the one that close to resonance with the atomic transition, e.g. 
the one with m = 1 with a maximum field at the lattice plane, the other modes hav- 
ing higher energies. In this part we consider isotropic optical lattice, hence cavity 
photons with a fixed polarization, i.e. TE and TM modes or any coherence superpo- 
sition of them. In Subsection 4.3 we treat an anisotropic optical lattice of uni-axial 
symmetry that give rise to photon polarization mixing. 

For the coupled exciton-photon system, in the dipole interaction and for the rotating 
wave approximation, the Hamiltonian is written as 

H = Yj n {^(k) B l B * + w ph( k ) a W + fk B l a k + ft a l B *\ • (26) 

k 

Here a^, a k are the creation and annihilation operators of a cavity photon with 
in-plane wave vector k. The cavity photon dispersion reads 



O) 



(27) 



where L is the distance between the cavity mirrors, and e ~ 1 is the cavity medium 
dielectric constant. The coupling parameter, e.g. for a lattice of square symmetry, 
is given by 



hfk = -i 



\tico ph (k)id 2 



(28) 



2e Q La 2 

In the case of excitons and photons close to resonance having small wave numbers, 
that is ti> ex (k) « co p h(k) for k ~ 0, the excitons can taken to be dispersion-less with 
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w e x(k) ~ u>a + 47. Translational symmetry leads to coupling between excitons and 
photons with the same in-plane wave vectors, namely to a momentum conservation. 

In the strong coupling regime where the exciton and photon damping rates are 
smaller than the coupling parameter, the Hamiltonian can be diagonalize by intro- 
ducing polaritons, through the transformation 

A k± = Xf B k + at, (29) 

where the exciton and photon amplitudes are 

(30) 

with Ak = ^5 2 k + \fk\ 2 , and we defined the exciton-photon detuning by 5k = (to p h(k)- 
a> ex )/2. In this regime the real system excitations are polaritons, quasi-particles 
formed by a coherent superposition of excitons and photons, with the system co- 
herently oscillating between those two. The polariton Hamiltonian is given by 

H = Y j nn r (k)Al r A kr , (3i) 

k; 

with the polariton dispersions 

OJph(k) + CO ex 

fl±(fc) = ^\ ± A.. (32) 

We get two polariton branches with a splitting of 2A^, the intersection point appear- 
ing at 5k = and the vacuum Rabi splitting equal to 2f k . The polaritons are denoted 
as upper and lower branches. 

The plots of Fig. 17(a) report the results for the upper and lower polariton branches 
in addition to the cavity photon and exciton dispersions, using the typical param- 
eters listed in the following. The exciton energy around small in-plane wave vec- 
tors is taken to be huj ex (0) = 2 eV, and a zero detuning between the exciton and 
the cavity photon dispersions at zero in-plane wave vector. The transition dipole 
is fi = 2 eh, and the optical lattice constant is a = 2000 A. The exciton-photon 
coupling parameter is |/| ~ 4 x 10 11 Hz, by neglecting the k dependence for small 
in-plane wave vectors. At zero in-plane wave vector the splitting is clearly given 
by the Rabi frequency. Fig. 17(b) reports the \X^\ 2 and \Y k \ 2 exciton and photon 
weights for the lower polariton branch. At the exciton-photon intersection point the 
polaritons are half exciton and half photon, that is \Xf\ 2 = \Yf\ 2 = 1/2. At large 
wave vectors the lower branch becomes excitonic, that is \X^\ 2 « 1, \Y^\ 2 « 0, and 
the upper branch becomes photonic, that is \X+\ 2 « 0, \Y£\ 2 « 1. 

The above discussion considered an ideal cavity with perfect mirrors. But in order 
to observe the system eigenmodes, we need to couple the internal cavity modes to 
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Fig. 17. (a) The upper and the lower polariton branches vs. the in-plane wave vector k. The 
dashed-dotted line denotes for the exciton dispersion, and the dashed-dotted parabola de- 
notes for the cavity photon dispersion, (b) The exciton and photon weights vs. the in-plane 
wave vector k for the lower polariton branch. Figure reprinted with permission from Zoubi 
and Ritsch (2007) Copyright 2007 by American Physical Society. 

the external world. This is done using the standard input-output formalism, which 
is well proved for high-Q cavities (Gardiner and Zoller, 2000). In this approach 
the internal cavity field modes and the external free radiation field are quantized 
separately and weakly coupled via the cavity mirrors. This leads to the cavity field 
damping but also allows to calculate consistently the relation between the cavity 
input and output fields. To include the atomic spontaneous emission as an energy 
loss term, we add a phenomenological effective exciton damping. This approach 
finally allows to calculate the transmission T, reflection R, and absorption A spectra 
for a given incident field in the cavity configuration of Fig. 16(a). 

For the previously investigated system we plot in Fig. 18 the transmission spectra 
for different wave vectors. Here we choose the following parameters for the damp- 
ing rates: exciton damping rate Y ex = 1 .5 x 10 9 Hz, upper and lower mirror damping 
rates y v = y L = 7.5 x 10 10 Hz- The peaks of the transmission correspond to the 
two polariton branches, which are the real eigenmodes of the system. The plot in- 
dicates that high transmission peaks are obtained for polaritons which are much 
more photonic than excitonic. At zero wave vector, as the polariton is half exciton 
and half photon, the upper and the lower spectra are identical. As the wave vector 
increase the upper branch becomes more photonic than excitonic, and the opposite 
for the lower one. Hence, the line-width of the upper branch becomes wider, which 
equals the cavity line-width, and the lower one becomes narrower, which equals 
the atom line-width. The reflection and absorption spectra can be easily derived in 
the same manner. The present results for the Mott insulator with one atom per site 
are expected to be significantly different from the other cases with other number of 
atoms per site and also from the superfluid phase. 

The widely used time-of-fiight method for observing the different quantum phases 
in optical lattices with their quantum phase transitions is a destructive observation 
tool (Greiner et al., 2002; Bloch et al., 2008). As the laser creating the lattice is 
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Fig. 18. The transmission spectrum vs. the angular frequency co for different in-plane wave 
vectors. The two peaks correspond to the upper and lower polariton branches at each in- 
plane wave vector. Figure reprinted with permission from Zoubi and Ritsch (2007) Copy- 
right 2007 by American Physical Society. 

switched off and after waiting some time a picture is taken, the system is destroyed 
and in each observation one need to prepare a new system. We suggest to use the 
optical linear spectra as a non-destructive observation tool for the quantum phases. 
The appearance of exciton and polariton provides a direct observation tool that can 
be repeated without modifying the quantum phase. Later we present other optical 
spectra, and the comparison between them allows us to extract information con- 
cerning the quantum phases. 



4.2.1 Two-Atoms per site Optical Lattice in a Planar Cavity 

The case of a two dimensional optical lattice in the Mott insulator phase with two 
atoms per site was investigated in Subsection 3.5. Here this system is introduced 
into a cavity where the optical lattice is taken to be localized in the middle and 
parallel to two planar cavity mirrors (Zoubi and Ritsch, 2008b), as in Fig. 16(a). 
The main objective is to examine how the previously derived dark and bright exci- 
tons behave within a cavity in the strong coupling regime. We start by considering 
the case of asymmetric optical lattice in which the two on-site atoms are sepa- 
rated by an average distance. The electronic transitions are coupled to the cavity 
photons with the interaction Hamiltonian Hi = £k,;,a Bfa k + fg*. <4-Bf ), wmi 
(or = 1,2) for the two atoms. The coupling parameter f" between atom a at site i 

and a photon with wave vector k is given by = -i yJfta> ph (k)ii 2 (2LNa 2 eb) _1 e' k I ? , 
where is the atom position. In terms of on-site symmetric and antisymmetric 

operators we get Hj = Zk,i,v (•/£■,• ^T a ^ + fk* a k^)' w ^ tn ( v = f° r tne two 
modes, where we defined the symmetric and antisymmetric coupling parameters 

fin = Wis + /k,i) and fk,i = (fit ~ fli> As tne on-site interatomic dis- 
tance is much smaller than the atomic transition wavelength, the symmetric state 
coupling is much larger than the antisymmetric one. The coupling can be controlled 
through the external laser fields creating the optical lattice in using different lasers 
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at different directions. It is possible to switch the antisymmetric coupling on and 
off by converting the lattice site from asymmetric into symmetric one. 

For symmetric optical lattices we have approximately /A = f£ t = f k j, where 

/ k s \ = V2/k,! and fg. = 0. Then the antisymmetric excitation states are almost 
decoupled from the cavity-photons, and they can be considered as dark states. Only 
the symmetric excitation states are strongly coupled to the cavity-photons. In this 
case by using the symmetric and antisymmetric excitonic picture in the momentum 
space, we get the total Hamiltonian 

H = Y J Ji {<(£) B^Bl + co ph {k) a{a k f k <a u + f k a^j + £ fc< BtfBf, (33) 

k i 



where f£ = f k = -i J £^ , and f£ = 0. In the strong coupling regime the 
symmetric excitons and the cavity photons are coherently mixed to form the new 
system excitations which are the polaritons. The above Hamiltonian first part can 
be easily diagonalize to give the following polariton Hamiltonian: 

H = ^ A[A kr + ^ BfBf. (34) 

kr i 

Here we can use exactly the definition of polaritons by replacing the exciton energy 
with the symmetric one and the coupling with the symmetric one. The antisym- 
metric states are on-site localized and dark, but they can be written and read in a 
controllable way, this fact making them useful as memory devices. 

4.3 Optical Lattices with Oriented Dipoles 

In the present Subsection we assume the transition dipole to have a fixed direction 
(Zoubi and Ritsch, 2009b). This case can occur for optically pumped atoms in 
a fixed Zeeman sub-level of a given angular momentum. Also molecules can have 
oriented transition dipole, e.g. diatomic molecules, where the electronically excited 
states have transition dipoles oriented along the molecule axis. In combination with 
external traps and fixed optical lattice laser polarizations the atomic and molecular 
transition dipole are aligned along a specific direction.This orientation that makes 
the lattice to be anisotropic with uniaxial symmetry, as represented in Fig 16(b) 
for a planar optical lattice. The anisotropic properties can be implemented as a 
signature for the formation of on-site molecules in the optical lattice. For example 
the formation of homonuclear diatomic molecules out of the two on-site atoms in 
the Mott insulator phase can be manifested through the anisotropic optical spectra. 
We will show that the electronic excitations and the cavity photons of both TE and 
TM polarizations are coherently mixed to form polariton branches, which can be 
directly observed in the linear optical spectra. For example, for an incident field of 
TE polarization we obtain both TE and TM transmitted and reflected fields. 
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The planar optical lattice is taken to be localized in the middle and parallel to cavity 
mirrors, as presented in Fig. 16(a). Here we generalize the previous presentation of 
the cavity photons in order to include the polarizations. The electromagnetic field is 
free in the cavity plane with in-plane wave- vector k, and is quantized in the perpen- 
dicular, z, direction with wave numbers k z = v j-, where L is the distance between 
the cavity mirrors, and m takes integer numbers, (m = 0, 1, 2, 3, • • • )■ Thus for each 
cavity mode (k, m) we have two possible polarizations TE and TM, which will be 
denoted by the indexes (s) and (p), respectively. The two cavity photon polariza- 
tions, (v = s, p), are assumed to be degenerate. Because we will consider only cav- 
ity modes with m = 1 to be close to resonance to the atomic transition, we have for 
the two polarizations o>k, 4 = &>k,p = co P h(k), with the cavity photon dispersion co p h(k) 
given by Eq. (28). The cavity photon Hamiltonian is H ph = Yik,vh°> p h(k) a kv a k , v , 
where and a kv are the creation and annihilation operators of a cavity photon 
in the mode (kv), respectively. The electric field operators and the field spatial 
functions are reported by (Haroche, 1992). We introduce the following unit vec- 
tors e z along the z axis, e k along the k in-plane wave vector, that is e k = k/fc, and 
% = 4 x e z , as illustrated in Fig. 16(b). 

The optical lattice is located at z = 0, the cavity middle, and parallel to the cavity 
mirrors. The coupling between the atomic transition and the cavity modes is given 
by the Hamiltonian H t = % J^ k v lf£ B k a^ v + f£* a^B^}, where the coupling param- 
eter is Hf£ = i -\Jj^- (m ' u v(k)). For example, we take jl to be real, e.g. for jl = [ix, 
see Fig. 21(a), and also we have e k = cos 6 x + sin 6 y and h k = sin 9 x, - cosQ y, to 
get %f k = iCu. sin 6 and %f£ = Ck (^) cos 6, where Ck = and oj = cn/L. 

The total Hamiltonian of the coupled excitations and photons is given by 




In the strong coupling regime the real system eigenmodes are cavity polaritons, 
which are obtained by diagonalizing the above Hamiltonian. We obtain the polari- 
ton Hamiltonian H = 2k, r ^1 A^-, characterized by three polariton branches, 

with the eigenfrequencies Q k± = " k ~l^" +A k and £lko = iOk, where A k = + \fk\ 2 

and |/ k | 2 = Yuv l/kl 2 ' w * m me ex citation-photon detuning 6 k = (oj k - oj ex )/2. For the 
case of exciton frequency of co ex /2n = 2.5 x 10 14 Hz, with a distance between the 
cavity mirrors of L = C7r/u) « 3.77 jim, zero detuning between the atomic tran- 
sition and the cavity photon at k = 0, at the angle 6 = n/4, transition dipole of 
p. = 2 eh, and lattice constant a = 2000 A = 0.2 ftm, Fig. 19 reports the three 
polariton frequency dispersions £l kr l2n as a function of k. 

In the limit of large detuning, where 6 k » |/ k |, we get Q. k+ ~ ^k + jg~, &k- ~ 
oi ex - and Q k o = w k . The upper (0) branch is a photon with a shift due to 
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Fig. 19. The three polariton frequency dispersions Q. r /2n vs. k, at 6 = n/4 and other param- 
eters reported in the text. Figure reprinted with permission from Zoubi and Ritsch (2009b) 
Copyright 2009 by Institute of Physics. 



the coupling to the excitation, and the lower (+) one is an excitation with a light 
shift due to the coupling to the cavity photon. The middle polariton branch is pure 
photonic. In that limit we get birefringence (Born and Wolf, 1997), determined by 
two refracted cavity fields, the ordinary field of the (0) branch, and the extraor- 
dinary field of the (+) branch. The polariton operators are in general a coherent 
superposition of excitations and photons of both polarizations, namely we have 
A kr = X kr B k + £y Y v kr a kv , with the relation \X kr \ 2 + £ v {Y^} 2 = 1, where the exci- 
tation and photon amplitudes of the upper (+) and the lower (-) polariton branches 



are X k± 



/ A k +. 
V 2A, 



2A r and Y L 



and the excitation and photon ampli- 



fi 



V2A k (A k +5 k )' 

tudes of the middle (0) polariton branch are X kQ = and Y£ = jh. For each photon 
polarization, in the strong coupling regime the longitudinal component along the 
transition dipole coupled to the excitation contributes to the upper and lower po- 
lariton branches. Instead the transverse component orthogonal to the dipole and 
decoupled from the excitation, creates the photonic middle branch. For each k, a 
direction exists where the polarization direction of a photon is orthogonal to the 
transition dipole, and hence no interaction between this photon and the material is 
obtained. 



To observe the system eigenmodes we need to couple the internal system to the ex- 
ternal world. Our observation tool is linear optical spectra. In assuming non-perfect 
mirrors, the internal system get coupled to their environment. For an incident ex- 
ternal field with a fixed in-plane wave vector and polarization we calculate the 
transmission, reflection, and absorption spectra, as shown in Fig. 16(a). The life- 
time of the excited state is included phenomenologically. For a TE or (s) polarized 
incident field, Fig. 20 report the transmission, reflection and absorption spectra of 
the (s) polarized fields, 7j , R s , and A (s \ respectively, as a function of frequency, 
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Fig. 20. Optical spectra vs the frequency in Hz for different angles, 6, at k = 5 x 10~ 7 A . 
In (a) the (s) polarized field transmission spectra, Tjf vs, of (s) polarized incident field; 
in (b) the (s) polarized field reflection spectra, Ry; in (c) the absorption spectra, A (s) ; in 
(d) the (p) polarized field transmission and reflection spectra, Tp and Rp \ Figure adapted 
with permission from Zoubi and Ritsch (2009b) Copyright 2009 by Institute of Physics. 



co — > co/2n, at different angles 9 and k = 5 x 10 7 A -1 . For the excitation and photon 
damping rates we used y/2n = 10 9 Hz and T ex /2n = 10 8 Hz. For 7^ and Rf the 
three peaks and three dips correspond to the three polariton branches. At 9 = we 
get zero transmission and complete reflection. The largest transmission peaks and 
the deepest reflection dips are obtained at = n/2. For the absorption spectra A (l) 
only two peaks are obtained correspond to the upper and lower branches, as the 
middle branch is pure photonic and no absorption take place, where the absorption 
is only for the polariton excitation part. Also here zero absorption at 9 = 0, and 
maximum absorption at 9 = n/2, are obtained. Fig. 20(d) reports the transmission 
and reflection spectra of the (p) polarized fields, Tp and R p s \ which are equal for 
identical cavity mirrors. Even though the incident field is (s) polarized we get trans- 
mission and reflection of (p) polarized fields, with maximum at 9 = n/4, and zeros 
at 9 = and 9 = n/2. 
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4.4 Finite Atomic Chain in an Optical Cavity 



In Subsection 3.2 we studied collective electronic excitations in a finite linear 
atomic chain, and in Subsection 3.4 we derived their damping rate into free space 
radiation. Now we localize the finite chain in the middle between two identical 
spherical cavity mirrors such that the chain axis is perpendicular to the cavity axis 
(Zoubi and Ritsch, 2009c), as shown in Fig. 14. Only the lowest Gaussian mode in 
the cavity is assumed to be close to resonance to the atomic transition. The electric 
field operator across the cavity waist, chosen along the optical lattice axis at z = 
and as a function of the distance from the cavity axis r, is given by (Haroche and 
Raimond, 2006) 

E(r) = i e~ r2/ < (e a - e* a*) , (36) 

where e is the photon polarization unit vector, vv the beam waist, and V the mode 
volume given by V = ttWqLo/4. Here E c is the cavity photon energy, and L is the 
distance between the cavity mirrors. 

The collective states are divided into two groups of symmetric and antisymmetric 
modes. The antisymmetric modes are found to be dark and decoupled from the 
cavity photons, while the symmetric ones are bright. For odd k-s, in the dipole 
interaction approximation, the coupling parameter is given by 



A = -M E ^ cot(^— ). (37) 

The lowest mode, k = 1, the one without nodes, has the highest coupling pa- 
rameter, that is nine times larger than the next bright mode of A: = 3. This mode 
dominates the optical properties of the system and can easily achieve the strong 
coupling regime within a cavity. It includes 81% of the oscillator strength, and is 
termed superradiant mode. Hence, in the following we consider only the coupling 
of such superradiant mode to the cavity photons. The Hamiltoian of the coupled 
superradiant mode and the cavity mode is given by the Jaynes-Cumming model 
H = E ex B^B + E c a) a + / B^a + f* a^B, where now the creation and annihilation 
operators B^ and B stand for the superradiant mode respectively, and a^ and a for 

the cavity photon. The coupling parameter is / = -i -yj ^(N+i) cot Gcyv+i) )' anc ^ me 
superradiant mode energy is E ex = E A + 2J cos ( j^f )■ The Hamiltonian can be eas- 
ily diagonalized using the transformation A+ = X ± B + Y ± a, to get H = ^E^ A \A r , 

with the two dressed eigenenergies E* = (E c + E ex )/2 ± A, where A = -^d 2 + \f\ 2 
and the detuning is 5 = (E c -E ex )/2. The eigenstates are coherent superpositions of 
the superradiant mode and the cavity photon. The superradiant mode amplitudes are 
X ± = + V(A + 6)/2A, and the cavity photon amplitudes are Y ± = // V2A(A + 8). 
At the intersection point, with 6 = 0, the dressed state is half a superradiant state 
and half s photon. But for large positive detuning the (-) branch tends to the su- 
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Fig. 21. The vacuum Rabi splitting frequency Qq vs. the atom number N, for interacting 
(full line), and non-interacting (dashed line), optical lattice ultracold atoms. The difference 
increases with the atom number. Figure reprinted with permission from Zoubi and Ritsch 
(2009c) Copyright 2009 by European Physical Society. 

perradiant state and the (+) branch becomes photon, and vice versa for the negative 
detuning. 

To emphasize the energy transfer influence within a cavity, we compare the re- 
sults to those without dipole-dipole interactions. The electronic excitation is de- 
scribed now by the Hamiltonian H ex = Yi n E A BlB n . The excitation-photon cou- 
pling Hamiltonian reads Hi = Ym (/ B\a + /* d*B n \ where the coupling parame- 



B = -j= Yin B n , to get the Hamiltonian H = E A B^B+E C a^a+f B^a+f* a^B, where 



/ = —i J ^f . In Fig. 21 we compare the vacuum Rabi splittings, £l Q = 2|/|//i, as 
a function of the atom number for the two cases, with and without the resonance 
dipole-dipole interaction. It appears that for a large atom number the vacuum Rabi 
splitting for independent atoms is larger than that of interacting atoms. We con- 
clude that the vacuum Rabi splitting is reduced by the dipole-dipole interactions 
and deviates from the square-root behavior. This result can be used as a signature 
for the formation of collective electronic excitations in atomic chains. 

4.5 One -Dimensional Optical Lattice coupled to a Tapered Nanofiber 

Recently new directions have been opened for the light-matter coupling in using 
optical fibers (Nayak et al., 2007). In tapered optical fibers the atoms are trapped 
outside the fiber and couple to the evanescent field surrounding the fiber. The guided 
modes of ultrathin optical fibers exhibit pronounced evanescent field that give rise 
to an array of optical microtraps. The system was realized recently for cesium atoms 



ter is / = -i J- 



. We define the collective excitation operator B n = -4= 
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Fig. 22. Two one-dimensional optical lattices with one atom per site, lattice constant a 
and distance d between them, are located parallel to a fiber, at the two opposite sides. An 
incident field is injected from the left side of the fiber and the transmitted field is observed 
on the right side. The transition dipole /u makes an angle 6 with the lattice direction. 

localized in one dimensional optical lattice parallel to the nanofiber (Vetsch et al., 
2010; Goban et al., 2012). We study here a simplified model related to this setup, 
based on a tapered optical nanofiber and two optical lattices (Zoubi and Ritsch, 
2010c). We consider two identical parallel one-dimensional optical lattices with 
one atom per site, separated by a distance d and located at the two opposite sides 
of the fiber, at a distance b from the fiber surface, as schematized in Fig. 22. 

The electronic excitation Hamiltonian in the two optical lattices is given by the 
one-dimensional version of the planar Hamiltonian of Eq. (18) 



Here (n,m) run over all the lattice sites, and (a, 8) stand for the two lattices (1,2), 
with B\ , B na are the creation and annihilation operators of an electronic excitation 
at site n in lattice a, respectively. The coupling parameter J"^ is for resonant dipole- 
dipole interactions, and gives rise to the excitation transfer among two lattice sites. 
The transfer term includes transfers among atoms in the same lattice, J l n ] n and J^f n , 
and transfers among atoms at different lattices, J X2 n and J 2X n . As before, the Hamilto- 
nian can be diagonalized in two steps. First, in the lattice sites in using the transfor- 
mation B na = -4= Yik e lkxL "Bka- Here x" is the position of site n in lattice a. Second, 
we diagonalize the Hamiltonian relative to the two lattice indexes, (or, 8), by ap- 
plying the transformation Bk v = Bk, ^ k2 . The new states represent entangled states 
between excitons from the two lattices, which are symmetric and antisymmetric 
states. The symmetric state is denoted by (v = s) and takes the plus sign (+), the an- 
tisymmetric state is denoted by (v = a) and takes the minus sign (-). The Hamilto- 
nian casts into the diagonal form H ex = v E v ex {k) B^ kv B kv , where the eigenenergies 
are E s ex (k) = E A + J(k) + J'(k) and E a ex {k) = E A + J(k) - J'{k). We define the exciton 
dynamical matrix by = £ L e ikL J a P(L), using fZ = J afi (L), with L = xi~ x%, 

where the dipole-dipole interaction is a function of the distance between the two 
atoms. For identical lattices, inside the same lattice we have J(k) = J n (k) = J 22 (k), 
and among different lattices we have J'(k) = J l2 (k) = J 2l (k). Excitons in one di- 
mensional optical lattices are presented in Subsection 3.2, and the coherent transfer 
parameter of excitons among two parallel lattices is discussed in Subsection 3.6. 




(38) 
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For the nanofiber photons we use a simplified picture. The one-dimensional optical 
fiber modes are given by the Hamiltonian H ph = Ytk^phik) The photon 

dispersion is taken to be of the form (o p h(k) = -^= ^Jk 2 + k 2 , where the wave number 

&o results of the fiber transverse confinement, and e is the fiber average dielectric 
constant. The electric field operator outside the fiber is 



E(r) = i J] -J-^T e u ( r) W ^ ~ a l e ~ ikZ ) ' < 39) 



where e is the photon linear polarization unit vector, V the normalization volume, 
and u(r) the mode function which includes the fiber field complexity. The exciton- 
photon coupling Hamiltonian is given by H in = J^ k {flf k ci^B, + %f£ a\B ks ], with 

the coupling parameter tif k = -iu{b) -J ^^f^ '■ Only the symmetric excitons are 
coupled to the fiber photons, where the antisymmetric ones are dark and decoupled 
from the photons. The total Hamiltonian reads 



H = ^ n [co a ex (k) B\B ka + (d'J® B\ B ks + co ph (k) a\a k 
k 

+ fk a k B\ s + f* a\B k ] , (40) 

where due to translational symmetry along the lattice and fiber axis, the Hamil- 
tonian is separated for each k. In the strong coupling regime we define the fiber 
polaritons by diagonalizing the Hamiltonian to get 

H = ^ fuo r pol (lc) A«Ai + J] *L*fa- (41) 

kr k 

We can use the polaritons definition of Subsection 4.2 in replacing the exciton and 
photon dispersions with the present one dimensional ones. 

Here we present the results for a system with the following parameters: lattice 
constant a = 1000 A; transition energy E A = 1 eV; transition dipole /u = 1 eh 
at an angle 6 = 90° with the lattice axis; distance between the two optical lat- 
tices d = 10a; fiber dielectric constant is e = 3; mode cross area S = Ana 2 , 
and mode function at the lattice position u(b) = 0.1. The cavity mode energy at 
k = is taken to be in resonance with the free atom transition energy, that is 
E p h(k = 0) = E A , where &o = E A ^fe/fic. The polariton dispersions are similar 
to those plotted in Figs. 17. The exciton-photon intersection point is here to at 
k = 0, where the polaritons are split by the Rabi splitting. To get the linear opti- 
cal spectra we consider an incident field from the far left side of the fiber, and we 
calculate the transmission spectra. The coupling of the fiber field to the external 
field at the two far edges of the fiber is included in the parameter y, and the fiber 
photon damping rate is taken to be Y ph which is assumed to be k independent. The 
symmetric exciton damping rate is included phenomenologically, and is taken to be 
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Fig. 23. (a) The transmission spectrum, for k - 10~ 6 A" 1 , and fiy = 10~ 4 eV. (b) A close 
look around the minimum. Figure reprinted with permission from Zoubi and Ritsch (2010c) 
Copyright 2010 by Institute of Physics. 



T s ex (k) as discussed in Subsection 3.4. The polariton damping rate is assumed to be 
r^ gl (k, 6) = T s ex (k) |X ± (fc)| 2 + T ph |y ± (fc)| 2 ; a small damping rate of fiT ph = 10" 10 eV 
is introduced for the photon, and the fiber edge coupling parameter is assumed to 
produce a bandwidth %y = 10" 4 eV. Fig. 23(a) reports the transmission spectra for 
a wavenumber k = 10" 6 A -1 , where at resonance we get a dip in the transmission 
spectrum. Fig. 23(b) representing a close look at the minimum shows a small shift 
relative to the free atom transition, produced by the dipole-dipole interaction. For 
this case the symmetric exciton line width is %Y s ex = 2.32 x 10" 8 eV, larger than the 
free atom line width of KT A = 2.5 x 10~ 9 eV. 



5 Optical Lattices with Defects: beyond the Mott Insulator State 



Up to this point we treated ideal optical lattices in the Mott insulator phase with a 
fixed number of atoms per site and in which all atoms are localized in the lowest 
Bloch band. In this section we partly relax this idealization and include different 
kinds of disorders into the system. In the following we introduce mainly two kinds 
of disorders. First, we exploit the excitation of on-site atoms into higher Bloch 
bands, where we concentrate in the lowest bands (Zoubi and Ritsch, 2009a, 2010a). 
Second, in the Mott insulator phase we consider localized defects that appear at 
several separated sites in the lattice (Zoubi and Ritsch, 2008a). We treat defects of 
vacancies distributed randomly in the lattice in the Mott insulator phase with one 
atom per site. 



40 



5. 1 Collective States of Atoms excited into higher Bloch Bands 



For deep optical lattices in the Mott insulator the lowest Bloch bands are well sep- 
arated and can be represented by vibrational modes of a quantum harmonic oscil- 
lator. The transition of atoms between different vibrational states can be induced, 
e.g., by inelastic Raman scattering of light, or by on-site atom-atom interactions and 
atom hopping among nearest neighbor sites. We present an alternative mechanism 
for the excitation and de-excitation of atoms between different Bloch bands through 
the coupling of the internal electronic excitation to the external atomic motion via 
electrostatic interactions (Zoubi and Ritsch, 2009a). We consider here the Mott in- 
sulator with one atom per site where the atoms are treated as two-level systems. 
In general ground and excited atoms in optical lattices experience different optical 
lattice potentials, and can have minima at different positions, the fact that induces 
the coupling of electronic excitations to the vibrational modes. We show that the 
on-site electronic excitation is dressed by a cloud of virtual vibrational modes that 
renormalize the excitation energy. In the appropriate regime, the excitation trans- 
fers among nearest neighbor sites due to resonance dipole-dipole interactions and 
can be accompanied by the emission or absorption of vibrational modes. 

We start from the electronic excitation Hamiltonian H = hiui^ - a>?) B]Bi + 
YiijhJij BjBj. The internal atomic transition frequency at site i of position n, is 
o)\ -o) 8 ., where we use the approximation cof = oj a + D^(uf" 1 ) for the ground and ex- 
cited state energies (A = e, g). Here Df n (uf m ) is the change in the local transition en- 
ergy due to an on-site atomic vibration excitation with (m > 0) vibrational quanta, 
and (jl> a is vibration-free transition frequency with zero quanta (m = 0). We take uf" 1 
to be the average spatial growth of the atomic size at site i in the A internal state due 
to the excitation of the atom to a higher vibrational state with m quanta. Df n (uf m ) 
is calculated relative to the local lattice ground state (m = 0), where D A (uf°) is 
included in o) x , and uf' n is measured relative to uf° of the vibration-less state. The 
excitation can exchange among atoms at different sites i and j due to dipole-dipole 
interactions, and it parametrized by the coupling parameter 7 ;; . This coupling de- 
pends on the average atomic vibration shift uf m and thus on the local vibrational 
excitation. In a deep lattice, uf m is a small deviation relative to the lattice constant 
a for the lowest vibrational modes. To lowest orders in this small perturbation we 
then can split the Hamiltonian into the form H = H ex +H vi b +H ex - vi b- Here H ex is the 
internal excitation Hamiltonian, which is obtained for atoms in the ground vibra- 
tional states, H vi b is the atom vibration Hamiltonian, for excited and ground state 
atoms, and the coupling Hamiltonian H ex - vi b is derived perturbatively for atoms 
excited to higher vibrational states. To the zero order, for atoms in the ground vi- 
brational states, we thus get H ex = ho) a B^Bi+Ytij B]Bp wnere the transition 
frequency o a = a)* - a>? is the same at each site, and TP. is the transfer parameter 
among atoms in the lowest vibrational states. This is exactly the Hamiltonian we 
treat before in Subsection 3.1, and that results in excitons for optical lattices. 
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Vibrational excitations of atoms in the ground and excited state potentials are ap- 
proximately described by the harmonic Hamiltonian 

H vib = £ fuJ v b]h + £ haj e v c\c u (42) 

i i 

where ool and of v are the vibration frequency for ground and excited state atoms, 
respectively, b], bi and cj, c ; are the creation and annihilation operators of a vi- 
brational mode at site i for ground and excited state atoms, respectively. The corre- 
sponding atomic displacement operators are 



H^ (i ' +i?M H^ (c ' +cI) ' <43) 

where m is the atomic mass. Here, ground and excited state atoms are considered 
as two different kinds of bosons, where each has its own optical lattice potential, 
as in Subsection 2.2. In the transition between the ground and the excited state one 
kind of bosons is destroyed and another created. 

To first order in the perturbation series with respect to uf, we get excitation- vibration 
coupling by the Hamiltonian H ex _ vib = H' ex _ vib + H' ex _ vib , where H' ex _ vib is for the on- 
site part, and H r J w for the transfer part. Processes in the first order terms include 
only a single vibrational quanta. To treat processes of more than a single vibrational 
quanta one needs to consider higher order terms of the perturbation series. As we 
treat here only the first order terms we drop the index (m) off Df n (uf m ). The on-site 
part is given by 

HL vib = X » [Mt [B] (B iCi ) + (Bjc]) B] - Mf B t + B] (*#)]} , (44) 

i 

where the coupling parameter is 



M -'fell^ • (45) 



H \dD A (uf) 

related to the slope of D A (uf). As D A (uf) is a function of the atomic wave function 
at site i and in a fixed Bloch band, the parameter Mf is related to the correlations 
between atomic wave functions at the initial and the final Bloch bands. The four 
possible transitions are presented in Fig. 24. The transfer part is given by 

H' e L, ih = Yj l \ F U c \ + F u b i + F u c J + F u b )] B < B J> (46) 
where the coupling parameter is 
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Fig. 24. The four possible transitions between the ground and excited electronic states 
including only single ground and excited vibrational states. Transition I includes a single 
excited vibrational quantum, and transition II includes a single ground vibrational quantum. 
Transition III is without any vibrational quanta, and transition IV includes ground and 
excited vibrational quanta. 

which is related to the derivative of Jy. 

In the limit of on-site excitation- vibration coupling which is stronger than the 
transfer coupling, namely in the limit of Mf » Ff 1 ., we can apply the following 

canonical transformation (Mahan, 1990) of O = e cr Oe~ cr and & = s B]B h with 

f ( b l - b ) - 5 ( c i - c ')\ to § et hi = bi ~% B < Bi and di = Ci + 5 B t Bi > 

and with B { = BjX , X = e~ s . The new effective local excitation represents an elec- 
tronic excitation dressed by a cloud of on-site virtual vibrations and can be con- 
sidered as an excitation-polaron. Namely, any electronic transition is accompanied 
by excitation and annihilation of ground and excited state vibrational modes. Here 
we assume the excitation-vibration coupling parameters to be site independent, by 
defining M A = Mf and Ff. = Ff 1 -. In terms of the new operators, the excitation 

Hamiltonian reads H ex = 2/^0*0 B^B, + J^jjtiJ® B]Bj, where coq = co a - A, with 

A = + ~~r~ • In this limit the effect of the on-site excitation- vibration cou- 
pling is a simple renormalization of the electronic excitation energy by a shift of 
A, where to a » A. This shift might be quite small but it needs to be considered 
for precision measurements as in a lattice clock transition energy (Takamoto et al., 
2005). The vibration Hamiltonian is now given by H vih = ficol b\bi + fioj e v etc;. 
The excitation- vibration coupling Hamiltonian of the transfer part is written as 
H ex -vib = Zij n [ F ij (ft + Cj) + F fj (h + SJ)] B]Bj. This term describes vibration 
changes induced by excitation transfer among different sites. Four possible pro- 
cesses are plotted in Fig. 25. They show the emission and absorption of a vibration 
due to the excitation transfer among nearest neighbor sites . 

In terms of transformed operators the interaction Hamiltonian can then be trans- 
formed into a momentum space representation in terms of exciton operators, with 
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Fig. 25. (a) Process I (full-line) for Ff. b]B; b] with an energy transfer accompanied by the 

ij i j j 

emission of ground state vibration at site j. Process II (dashed-line) for Ff. B;Bj c] with an 
energy transfer accompanied by the emission of excited state vibration at site i. (b) Process 
III (full-line) for Ff fiffi, b{ with an energy transfer accompanied by the absorption of 

ground state vibration at site i. Process IV (dashed-line) for Ff. B^Bj cj with an energy 
transfer accompanied by the absorption of excited state vibration at site j. 

the exciton Hamiltonian H ex = £ k ha>(k) B^B k , for two dimensional optical lattices. 
The vibration operators, even though representing on-site localized vibrations, can 
be rewritten equally in the momentum space by applying b t = -4= 2 q e iqn 'bq and 

Ci = -jfi Zq e iqDi c q , to get H vib = £ q fio/ v blb q + 2 q hu e v c\c q , which have flat dis- 
persions, namely ^-independent. The exciton- vibration coupling now reads 



H ex - Vib =J]h [F e (k + q) c q + F*(k) b q + F e (k) c! q 

k.q 

+ FZ(k + q)E q } Bl +q B k , (48) 

where F A (k) = -j= F A (L)e' k ' L . The Hamiltonian H ex - vi b describes scattering of 
excitons between different wave vectors involving the emission and absorption of 
a ground or excited state vibrational quantum. Such a process is possible in the 
limit of J » to A , where the exciton band width is larger than the vibration energy. 
The scattering conserves energy and momentum, and as the vibration dispersion 
is a flat one, the vibration can absorb any amount of momentum from the exciton. 
The exciton transition rate within a spontaneous emission is wj* oc \hF A (k)\ 2 , which 
results of the Fermi golden rule. The exciton- vibration interaction will serve as an 
important source for relaxation of excitons toward lower energies of the exciton 
band and therefore a thermal equilibrium in optical lattices. 

The exciton- vibration coupling yields polariton- vibration coupling (Zoubi and Ritsch, 
2010a). In using the inverse previous transformation B k = J^ r Xl* A kr , we get the 
polariton- vibration interaction by 
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H pol - vib = ^ 2j n \ Fe( - k + q) C q + ^ K + ^( k ) C -q 

+ F*(k + q) *y (xr^ +?) ) A[ k+q) A kr , (49) 

where the exciton amplitudes indicate that the interaction is due to the polariton 
excitonic parts. This interaction describes scattering of polaritons between states 
with different wavevectors by the emission or absorption of vibrational quanta. As 
the vibrations are dispersionless, one need to take care only of the energy conserva- 
tion. A scattering takes place between different momentum states in the lower and 
the upper polariton branches. But if the vibration energy equals the splitting energy 
between the two branches, the polaritons can jump between the two branches by 
the emission or absorption of vibrations. The polariton- vibration interaction can 
represent a significant source for polaritons relaxation toward the minimum energy 
at k = 0. Using the Fermi golden rule, the damping rate for the spontaneous emis- 
sion of vibrations at zero temperature T = off polaritons at the lower branch is 

given by oc (^(k)! 2 + \F s (k - q)| 2 ) |X~| 2 X~ k _ q) . The appearance of the ex- 
citon amplitudes indicate that the process is much more efficient in regions where 
the polaritons are more excitonic than photonic. 



5.2 Excitons and Polaritons scattering by Defects in the Mott Insulator 



The appearance of some defects is unavoidable due to imperfections in the dynam- 
ical formation of the Mott insulator phase, in particular there can be missing or 
extra atoms at some sites. We investigate here the effect of such defects on the dy- 
namical properties of excitons and cavity polaritons (Zoubi and Ritsch, 2008a). We 
concentrate on the case of a low defect number, where the exciton and cavity polari- 
ton picture still holds and these quasi-particles are only scattered by single defects. 
Namely, the distance between each two defects is large enough for the formation 
of coherent excitons and cavity polaritons, which propagate as free quasi-particles 
between each two scattering processes. In principle the defects can move by hop- 
ping among the lattice sites, but the corresponding time scale is so long that they 
can be considered frozen. Here we concentrate on one type of defect which is in- 
duced by a missing atom, that in the case of one atom per site corresponds to the 
presence of vacant sites. In calculating the exciton and cavity polariton elastic scat- 
tering amplitude of such defects, we show that vacancies behave like hard sphere 
scatterers. We suggest cavity polaritons as a tool to detect such defects and show 
how the scattering can be controlled by changing the exciton -photon detuning. 

We treat a two dimensional optical lattice filled with one atom per site in the Mott- 
insulator phase. A single atom missing at the origin (r, = 0) then creates an impurity 
in the artificial lattice of ultracold atoms, which for sufficient lattice depth stays lo- 
calized and will not hop among the lattice sites as shown in Fig. 26. We study the 
scattering of excitons off such a missing atom. The addition of the impurity to the 
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Fig. 26. Schematization of a vacancy in one dimensional optical lattice. 

ideal atom crystal needs only a small change to the system Hamiltonian H = H + V 
presented in Subsection 3.3, whose components now read H = Y,^fioj ex (k) 5^5 k 
and V = -tiu) A bLBq, where H represents free excitons, and u> ex (k) is the exci- 
ton dispersion with in-plane wave vector k. V is the impurity Hamiltonian at the 
origin. To get the above form, we add and subtract the impurity Hamiltonian V to 
the whole Hamiltonian H, that is H = H - V + V, and then we define the ideal 
case Hamiltonian by H = H - V. For an exciton the impurity thus appears as a 
potential well located at the origin, with depth %to A and radius a. We will not try to 
find the self-consistent eigenstates of H and consider only the scattering problem. 
Because we can neglect trapping of an exciton in the impurity potential as a vacant 
site cannot absorb the trapping energy, we consider only the scattering process of 
an exciton off the impurity and calculate its scattering amplitude. The initial exciton 
initially with wavevector k very far from the impurity, is scattered elastically into 
the wavevector k' also very far, with |k'| = |k|. As the impurity is very deep with 
depth E A = h<x>A, in order to calculate the scattering amplitude we need to proceed 
beyond the Born approximation even though the perturbation is confined to a ra- 
dius a. We apply the Schwinger-Lippmann equation (Newton, 1982) in the case of 
small wave vector scattering excitons, that is ka «: 1. For this case we can use the 
parabolic approximation for the dispersion of the excitons, co ex (k) = u) ex (0) + , 
with m ex the exciton effective mass, and for isotropic atoms in a square lattice of 
cubic symmetry as in Subsection 3.3 m ex = -h/(2Ja 2 ). We obtain the scattering 
amplitude 

f(k) = 7^i—, t, (50) 

i + £[*>(*)-'!] 

where the effective exciton band width is defined by A ex = ^f al ■ The dipole- 
dipole interaction energy between different sites in the optical lattice is small, and 
hence the excitons have a small band width. In the present case we have the limit 
of E A » A ex so that the scattering amplitude is f(k) « fa/fe) > wn i cn exactly repro- 
duces the result for the scattering off a hard disk of radius a with the scattering cross 
section defined by cr(k) = 2n\f(k)\ 2 . We thus conclude that an impurity generated 
by a missing atom in an optical lattice acts effectively like a hard disk of radius a. 
Hence if a large number of incident excitons with identical wave vector are scat- 
tered elastically off the impurity, we get a ring of radius k of scattered excitons. The 
scattering amplitude is plotted in Fig. 27 as a function of wave vectors, k, for lattice 
constant a = 2000 A. The singularity at k = is the two dimensional signature. 
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Fig. 27. The scattering amplitude vs. wave vector k, for zero detuning and a lattice constant 
a = 2000 A. Figure reprinted with permission from Zoubi and Ritsch (2008a) Copyright 
2008 by Institute of Physics. 

We now add cavity mirrors to our lattice as described in Subsection 4.2, where in 
the strong coupling regime, the system eigenstates are the cavity polaritons. As a 
polariton is a coherent superposition of an exciton and a photon it can be scattered 
off an impurity due its excitonic part. The corresponding Hamiltonian then reads 
Ho = Jj kr HQ r (k) aLA^t and V = -E A B^Bo, where the upper and lower polariton 
branches are presented explicitly in Subsection 4.2. We calculate the polariton scat- 
tering off such impurity. Again the scattering is elastic so that an incident polariton 
in branch r with wave vector k will scatter into a polariton with wave vector k' in 
the same branch. We consider the scattering of small wave vector polaritons, that 
is in the limit ka «; 1 in the lower branch. In this limit the polariton dispersion can 
considered approximately parabolic with a polariton effective mass of m p . This is 
of the order of the cavity photon effective mass m p « (Hn)/(cL). Hence, the lower 
polariton branch dispersion is taken to be E-{k) = £-(0) + j*-, and we have in 
general m p «; m ex . The scattered states will also have small wavevectors due to 
the energy conservation. As the upper branch has higher energies its contributions 
are negligibly small. In addition the upper branch for larger wavevectors is mainly 
photonic with small excitonic part, and thus it contributes only weakly to the im- 
purity scattering. We consider the case around zero detuning between the excitons 
and photons. Therefore the excitonic weight \X^\ 2 will change from half around 
zero wave vector up to one for large wave vectors, where the lower branch became 
excitonic. Because we neglect the contribution of the upper polariton branch, we 
drop the branch index, and all the parameters will be only for the lower branch. For 
simplicity, we use a model for the lower branch dispersion in place of the real one. 
The lower branch is divided into two parts: the first part between < k < k 
is taken to be of a parabolic dispersion with a polariton effective mass, where 
E(k) = E(0) + (Ti 2 k 2 ) i '(2m p ); and the second part between k < k < n/a, where 
n/a is the Brillouin boundary, is taken to be dispersionless with energy equal to 
the exciton energy at zero wave vector, with k <S n/a. The intersection point be- 
tween the two parts is fixed by E = ti 2 ky{2m p ) + E(0). The calculation yields the 
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Fig. 28. The scattering amplitude vs. exciton-photon detuning, for k ~ polaritons. Figure 
reprinted with permission from Zoubi and Ritsch (2008a) Copyright 2008 by Institute of 
Physics. 

scattering amplitude 




where we defined the polariton effective band width A„ = and A* = E -E(k). 
As E A > A k , we have (nE A ) /(4A k ) > 1, and the scattering amplitude is negative, 
that is f(k) < 0. Thus the above impurity effective potential is repulsive for the 
polaritons. In the case of zero detuning, that is do = 0, for small wavevectors \X k \ 2 = 
1/2, Ak of the order of the exciton-photon coupling, \g k \, where E A » A*, we 
get f(k) ~ -A k /A p . This scattering amplitude equals that of the scattering from 
an effective potential of a square barrier potential of height A* and width a. For 
the case of negative detuning, in the limit of E A » So, where A k is of the order 
of the exciton-photon detuning, the scattering amplitude is f(k) « -(2X 2 A k )/ A p . 
This scattering amplitude equals that of the scattering from an effective potential 
of a square barrier potential of height 2X 2 A k and width a. Hence we find that the 
scattering amplitude can be controlled by changing the exciton-photon detuning. 

Using the following parameters: atomic energy E A = 2 eV, exciton-photon cou- 
pling fi\g\ = 0.0001 eV, lattice constant a = 2000 A, approximately normal incident 
waves with k = 10" 6 A, the scattering amplitude is calculated as a function of the 
exciton-photon detuning, as reported in Fig. 28. It is clear that maximum scattering 
is obtained for zero detuning. As the detuning increases, either positive or negative, 
the k « polaritons become more photonic and they are much less scattered off the 
impurity. As the atomic energy E A is much larger than the exciton-photon coupling 
energy %g k , the scattering amplitude shows only a negligibly small asymmetry for 
negative and positive detunings, but we expect a clear resonance of the scattering 
amplitude around a resonant value 8q. Considering a stream of incident polaritons 
with the same wave vector k = fckwe then should see a ring of radius k of scattered 
polaritons. The effect of defects on the dynamics of excitons in optical lattices can 
be generalized to the case of the Mott insulator phase with several atoms per site 
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that contain defects of extra or less atoms at specific sites, as we presented for the 
case of two atoms per site with a defect of a singly occupied site. 



6 Conclusions 

Excitons and polaritons in an optical lattice filled with ultracold atoms in the Mott 
insulator phase represent an idealized form of an optically excitable solid state crys- 
tal. Modeling by a two component Bose-Hubbard model allows us to derive the 
quantum phase diagram identifying the regions of Mott insulator phases. Excitons 
occur as natural elementary excitations in such a system including energy trans- 
fer induced by resonant dipole-dipole interaction. We extract the condition for the 
appearance of excitons in such artificial crystals by comparing between the trans- 
fer parameter and the excited atom decay. Excitons in the Mott insulator phase 
are collective delocalized electronic excitations in the lattice with their momentum 
distributed among many lattice sites. Propagating excitons are associated to large 
systems, while systems with boundaries require standing wave excitons. In general 
their damping rates are large and deviate from single atom decay, where excitons 
can be metastable or superradiant depending on their wave vectors and polariza- 
tions. Dark or metastable excitons open the door for storage and manipulation of 
photonic qubits with a wide range of applications for quantum information pro- 
cessing. Coherent transfer of excitons among different optical lattices differs from 
hopping of electronic excitations among two separated atoms. Here the transfer 
amplitudes decay exponentially with the inter-lattices distance. 

Lattice excitons can be used as active medium for a cavity QED in the resonant 
regime. Cavity photons affect the boundaries between quantum phases and in these 
systems the Mott insulator phases appear only for deeper optical lattices. This gives 
new tools to observe and control the quantum phase transitions via the average 
number of cavity photons. In the strong coupling regime the excitons and pho- 
tons are coherently coupled to form new quasi-particles termed polaritons. Optical 
linear spectra are presented as a nondestructive observation tool for the system 
different quantum phases. Systems of various dimensionality and geometries are 
investigated, ranging from a finite chain between spherical cavity mirrors up to a 
two-dimensional optical lattice between planar cavity mirrors. Moreover, the for- 
mation of molecules in optical lattices results in an anisotropic crystal that can be 
manifested through cavity a photon polarization mixing in the linear spectra. We 
presented also a set-up of two parallel one-dimensional optical lattices formed by an 
optical nanofiber and that provides a tool for a direct excitation of one-dimensional 
excitons through the coupling to fiber photons. Moreover the system gives rise to 
one-dimensional polaritons that dominate the optical properties in the strong cou- 
pling regime. We suggested a new mechanism for the excitation and de-excitation 
of atoms into higher Bloch bands through the energy transfer process, and we pre- 
sented it as a significant source for the thermalization of excitons and polaritons. 
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The elastic scattering of excitons and polaritons off vacancies are used for verifying 
the presence of various defects that can appear naturally in optical lattices. 

The discussion in the present review though presented for optical lattice ultracold 
atoms can be adapted directly for various systems containing lattices of active ma- 
terials. For example, one can imagine a lattice of semiconductor nano- structures, 
e.g., chains of quantum dots, arrays of color centers in solids, organic molecules 
ordered on a matrix, and Rydberg atoms in a lattice. The results can be extended 
into various directions, just to mention the implementation of excitons and polari- 
tons for the physical realization of quantum information processing. The parameter 
controllability and the wide range of lifetimes of excitons make them a strong can- 
didate as qubits. But efficient quantum processing implies entangled qubits, which 
are achievable via interaction among such qubits. Therefore our results should be 
extended to exploit exciton-exciton interactions through saturation and Coulomb 
forces, which also can find application for photonics using optical lattices. We 
expect this research to guide future experiments of electronically excited optical 
lattice ultracold atoms both in free apace and for cavity QED. 
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